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lift slope, or nondimens ional distance between midchord of the wing and 
elastic axis (positive aft) 

semichord, c/2 

coefficients in equations (73) and (80) 
lift deficiency function, equation (B6) 
airfoil profile drag coefficient 

horizontal force, side force, and thrust coefficients 

rolling, pitching, and yawing moment coefficients 

lag damper coefficient 

chord length 

induced drag 

expectation of ( ) 

distance between elastic axis and feathering axis at blade root 
distance between mass center and elastic axis of blade cross section 
aerodynamic loads on the wing section, equation (Bll) 
net hub forces equal to (H,Y,T) 
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hub forces vector of the £th blade, (F ,F ,F )„ 

x x z % 

gust frequency 

gravitational acceleration or the gust velocity vector 
horizontal force in the rotor frame = p (Rfi) 2 (ttR 2 )Cjj 
plunging motion of wing cross section 
moment 

imaginary unit = /-I 

shaft inclination angle (positive fore) 
quadratic performance function, equation (72) 

Bessel function of the nth order 
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of inertia of a blade about flapping hinge = J r (r - rg) m dr 
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feedback gain vectors, equation (26) 

polar radius of gyration of a blade per unit span 

Kalman gain at n time-cycle 

blade lift per unit span 

blade pitching moment per unit span 

aerodynamic moments on the wing section, equation (Bl2) 
variance of the error before the measurement at n time-cycle 
net rolling moment = p (Rfi) 1 (TrR 2 )RCj^ 
net pitching moment = p (Rfi) 2 (ttR^RC^y 
net yawing moment = p (Rfi) 2 (irR^RC^ 

J > R 

r(r - r fi )m dr 
r e V 

hub moments vector of the £th blade, (M x >My,M z )'^ 
mass per unit span 

generalized mass of the jth mode, equation (2a) 
number of blades 

angular velocity of a rotor along the Xg-axis 

variance of the error after the measurement at n time-cycle 

angular velocity of a rotor along the Yg-axis, Floquet transition 
matrix, variance of the process noise, or hub inertia moment vector 
in equation (B3) 

generalized force of the jth mode, equation (2c) 
generalized coordinate of the jth mode 

kth constant coefficient in the solution of Aq, equation (7) 
perturbation generalized coordinate vector, equation (5) 
rotor radius 

angular velocity of rotor along the Zg-axis 

radial coordinate of blade or variance of the measurement noise 
position of flapping hinge 
position of lag damper 
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thrust in the rotor frame = p (Rfi) 2 (ttR 2 ) C^, or period 
transfer function matrix 
matrix, equation (59) 
time 

velocity of rotor along the Xg-axis 

tangential, normal, and radial components of velocity at a blade element 

radial deflection of a blade; process noise 

horizontal component of gust velocity 

velocity of a rotor along the Yg-axis 

lead-lag deflection of a blade; measurement noise 

lateral component of gust velocity 

velocity of a rotor along the Zg-axis; gross weight 

weighting matrices, equation (72) 

flap deflection of a blade 

vertical component of gust velocity 

gust amplitude 

precone angle 

nondimens ional radial coordinate = r/R 

side force in the rotor frame = p (RH ) 2 (iTR ? ')Cy 

measurements 

angle of attack 

Lock number = pacR^/lg 

perturbation from a steady value, or small increment 
Kronecker delta function 

blade trim pitch input = 0^ + 0 fc (x - 0.75) + |0 cos tp + 0 lg sin if/ 
blade twist rate 
collective pitch control 
lateral pitch control 
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Subscripts : 


longitudinal pitch control 
control vector, equation (59) 
control inputs, equation (57) 
inflow ratio = (Vg + V sin ig)/Rft 
advance ratio = U cos ig/Rfi 
mean induced velocity 
air density 
solidity = Nc/ttR 

torsional deflection of a blade, inflow angle 

phase shift in figure 17 

azimuthal angle 

rotor rotational speed 

gust angular velocity = 27rfQ 

jth natural frequency of a blade 

lag frequency of a blade 


or phase angle 
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body frame 
blade frame 

control or cosine element 
gust 

rotational hub frame 
inertial frame 
jth mode 
&th blade 
n time-cycle 
rotor frame 
sine element 

initial value, amplitude, uncontrolled value, or origin 
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SUMMARY 


Two control schemes designed to alleviate gust-induced vibration are analytically 
investigated for a helicopter with four articulated blades. One is an individual 
blade pitch control scheme. The other is an adaptive blade pitch control algorithm 
based on linear optimal control theory. In both controllers, control inputs to 
alleviate gust response are superimposed on the conventional control inputs required 
to maintain the trim condition. A sinusoidal vertical gust model and a step gust 
model are used. The individual blade pitch control, in this research, is composed 
of sensors and a pitch control actuator for each blade. Each sensor can detect flap- 
wise (or lead-lag or torsionwise) deflection of the respective blade. The actuator 
controls the blade pitch angle for gust alleviation. Theoretical calculations to 
predict the performance of this feedback system have been conducted by means of the 
harmonic method. The adaptive blade pitch control system is composed of a set of 
measurements (oscillatory hub forces and moments) , an identification system using a 
Kalman filter, and a control system based on the minimization of the quadratic per- 
formance function. Calculations for the individual blade pitch control system show 
that the thrust fluctuation can be reduced by more than 50% in the low gust frequency 
range. In addition, it was found that the adaptive blade pitch control system can 
also be useful for high frequency gusts. 


1 . INTRODUCTION 


The helicopter has the capability of flying close to the ground where the motion 
of the atmosphere may be thought of as turbulent flow. To avoid structural vibra- 
tions and unfavorable dynamic characteristics in flying and riding qualities caused 
by gusty winds, it is necessary to analyze the gust response of a rotary wing air- 
craft and to develop a scheme to alleviate such responses. 

Studies on the gust response of a helicopter rotor are reported in refer- 
ences 1-9. Many different assumptions were considered in those analyses: for exam- 

ple, quasi-steady or unsteady aerodynamics; constant or variable induced velocity dis- 
tribution; rigid or flexible blades, with or without hub motion; and sudden or gradual 
penetration into the gust. In references 5 and 9 the wind tunnel response of a rotor 
encountering a vertical gust was experimentally investigated. Experimental data for 
vertical gusts show that the rotor thrust response is most sensitive to the gust 
amplitude and frequency. For the sinusoidal gust with frequency ojq, the thrust 
response characteristics have not only a frequency Wq component, but also 
(o)q + nft) (n = ±1, ±2, . . j.) components in the fluctuation, where ft is the rotor 
rotational speed. In this thrust response, the Wq component of the oscillatory 
characteristics is dominant in the comparatively low gust frequency range (coq << ft) . 
Other researchers have investigated the influence of atmospheric turbulence on rotor 
aeromechanical stability. In references 10-18, rotor stability analyses are reported 
for the case of a rotor encountering a gust. Typically the turbulence is a random 
gust with stochastic properties. 



Attempts to alleviate gust-induced vibration have been made in the work 
described in references 19-37. Among them, Briczinski and Cooper (ref. 19) and 
Briczinski (ref. 20) investigated the effect of a rotor/vehicle state feedback sys- 
tem on the handling qualities of a helicopter, specifically the characteristics con- 
cerned with gust response. It was found that the feedback scheme of the rotor tip- 
path-plane or body-state was very useful as a means of gust suppression. Frick and 
Johnson (ref. 21) and Johnson (ref. 22) studied the performance of an optimal control 
system applied to proprotor/wing response to vertical gusts. In that investigation, 
the von Kdrmdn model or the Markov-process gust model were assumed as gust models, 
and linearized state equations were assumed to govern the motion of the rotor and 
wing for the rotor model. Optimal feedback inputs were determined by using linear 
optimal control theory, and state variables were estimated by a Kalman-Bucy filter. 
Significant and simultaneous reduction in the rotor and wing response was achieved. 
Ham et al. (ref. 23) and other investigators (refs. 24-26) have dealt with a classic 
parameter optimization method to alleviate the gust response of the tilt rotor air- 
craft. A von Kdrmdn gust model was also used in these studies. Zwicke et al. 

(ref. 27) and Taylor et al. (ref. 28) investigated the performance of the optimal 
sampled-data feedback system on the gust response. Zwicke also studied a suboptimal 
feedback system derived from the above control system; a significant reduction in the 
gust response was achieved by using this suboptimal feedback system. Ham and 
McKillip (ref. 29) and Ham (ref. 30) developed an individual blade control (IBC) 
design for gust response alleviation, which has also been, applied to suppressing 
blade stall-induced vibration, lag damping augmentation, and other dynamic problems 
(refs. 31-36). Saito et al. (ref. 37) have also studied a simple feedback system to 
alleviate rotor gust response. In their control scheme, individual blade pitch angle 
control, using scheduled feedback gains derived from analytical calculation, was used. 

There has been significant progress in vibration reduction techniques for heli- 
copters, theoretically and experimentally, in the past decade (refs. 37-83). In 
these vibration reduction systems, control schemes, known as multicyclic (ref. 45), 
or higher harmonic control (ref. 61), were applied to reducing the inherent vibratory 
response of a helicopter. Pitch inputs at harmonics of the rotor rotational speed 
are used. Typically, the helicopter is represented by a linear, quasi-static fre- 
quency domain model. The relationship between control inputs, 0, and outputs, z 
(which can include loads, vibrations, and rotor performance parameters) are modeled 
by a transfer matrix Tq. Theoretical and experimental results show that the level 
of vibration in a helicopter can be significantly reduced by using controllers in 
which the transfer matrix is updated by a Kalman filtering scheme. These vibration 
reduction systems are reviewed in detail in reference 70. Investigations have 
recently been extended to the analysis of nonlinear effects on controller perfor- 
mance (refs. 71 and 72). Other researchers investigated structural modification of 
the rotor blade or the fuselage body for vibration reduction (refs. 79-83). 

In this report, a gust alleviation system is studied using two control schemes: 
first, an individual blade pitch control scheme, and second, an adaptive blade pitch 
control scheme based on linear optimal control theory. The emphasis is on the appli- 
cation of these controllers to alleviate gust response as measured by blade flapping 
and hub oscillatory forces and moments. In this study, the local momentum theory 
(LMT) (ref. 84) is used to calculate timewise vibratory airloadings and moments at 
the rotor hub. The LMT is based on the instantaneous balance between the fluid 
momentum and the blade elemental lift at a local station in the rotor rotational 
plane. Therefore, this theory can be used to evaluate timewise variation of airload- 
ing, aerodynamic moments, and the induced velocity distribution along a blade span. 

The LMT has been applied to study many rotary wing phenomena (refs. 9, 37, 85-87). 

An H-34 rotor model with four articulated blades is used. The blades have full 
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flexibility in the flap, lead-lag, and torsion directions. The perturbational 
motion of the fuselage is not considered, and the shaft inclination angle ig is 
assumed to be zero (perpendicular to the flight direction) . 

The author expresses his great gratitude to Dr. Wayne Johnson and Dr. William 
Warmbrodt for their review and valuable suggestions. Further thanks are extended to 
the late Mr. John L. McCloud III and to Mr. Stephen A. Jacklin for their kind and 
encouraging discussions. The author also acknowledges computer operation support by 
the people of NASA Ames Research Center. This study was made possible by the National 
Research Council Associateship . 

2. INDIVIDUAL BLADE PITCH CONTROL 


In this section, an individual blade pitch control scheme is studied as a gust 
alleviation system. First, the harmonic method is applied to the theoretical calcu- 
lations of the rotor gust response with a simple feedback control system for blade 
pitch. A similar feedback control system for an individual blade was investigated 
for a rigid blade in reference 37. The analysis is extended here to include elastic 
blades; therefore, sensors to detect blade deflection (flap, lead-lag, torsion) are 
considered to study the sensitivity of the control system. A blade stability analy- 
sis with the feedback control system was made using Floquet theory. Finally, for a 
specified feedback gain, this control system is validated for a rotor penetrating 
into a vertical gust. 


Theoretical Calculation by the Harmonic Method 


Theoretical calculations of rotor response t'o vertical gusts in forward flight 
have been made by using the harmonic method. Derivation of the isolated blade equa- 
tions of motion and method of solution are based on references 88-95 with the excep- 
tions that the hub motion (U, V, W, P, Q, R') is replaced by a vertical gust velocity 
w G . Coordinate systems used in this study are briefly explained in appendix A. The 
displacements of an elastic blade are expressed as follows: 


" - E Sjqj . 
J = 1 


V = ^ ^j q j ’ 
J=l 


<t> = £ 

j=i 


(i) 


where qj is the generalized coordinate of the jth coupled mode and w j , v j , and 
<j)j are the corresponding jth coupled mode shapes. When the associated jth modal 
frequency is denoted by toj , the following equation can be derived according to the 
Rayleigh-Ritz approach (ref. 96): 

m j q j + “j^j = Q j (j = 1 > 2 > 3 > * • '•) ( 2 a) 


where 


”, ■ r 


t(k <(>j + eWj - e0Vj)<j)j + (w^ + ecf>..)w.. + (vj - ee^v^jm dr 


(2b) 
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(2c) 


Qj = J ({M ax - m[eg + k 2 (0 + Q 2 Q + 2ft0v’) + e0(e o Q 2 - 20u)]}$j 

o 

+ [F az - m(g + e0 - eft 2 0) - 2mfiWQv]Wj + {F ay + m[-2f2u + 2eft(v* + 0w') 
+ e00 + fi 2 (e Q + 2e)] - 2mfiu}vj^dr 


In this study, the radial extension of the blade is neglected in the blade equations. 
However, to take account of the Coriolis force caused by the flapping, u is approxi- 
mated as 



r 


[1 + (vq + v') 2 + (wq + w’)^]dr^ - j r 




J (v' 2 + w ,2 )dr - J (v a v' + w a w’)dr 



(1 + v' 2 + w ' z ) dr 


| 2 > 


r 


00 

■ ? 

J = i 


(3a) 


where 



WqW' dr 


(3b) 


and Wq is the time-averaged flap displacement. The lead-lag displacement is 
omitted in equation (3). More detailed expressions of M ax , F az , F ay , etc., shown 
in equation (2), are given in appendix B. For this study, consider the vertical gust 
component Wq as an external exciting source. By introducing wq into the aero- 
dynamic loads and eliminating W in reference 89, the blade equation of motion, 
including gust excitation, can be obtained. Equation (2) yields a set of simulta- 
neous differential equations for q j . A trimmed value of qj can be obtained by 
setting all perturbing values from trim condition to be zero. When the gust velocity 
defined in the rotor frame 


W G “ w Go exp[i(co G t + <f>) ] (4) 

is superimposed on the trim condition, additional blade motion will occur, where w Go 
is a constant real number. If the gust-induced blade response amplitude is assumed 
to be small, the incremental blade motion can be described by a small perturbation 
Aqj from the trim value qj . Furthermore the perturbation generalized coordinate 
vector, Aq, is defined as 


Aq (Aq x , Aq£, Aq^, . • .) (3) 

It is then easy to show that Aq can be described by a set of second-order linear 
differential equations. Hereinafter, all quantities (i.e., elastic deflections, 
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aerodynamic forces, and moments, etc.) should be considered as perturbed values from 
a trim position. 

In forward flight, coefficients of the differential equations become time- 
varying and can be written as follows: 

3 

Z3 { [Ajjj exp(Jmi(;)]Aq + [B m exp (imifj) ] Aq + [C m exp (inu|>) ] Aq} = £ F n exp[i(w G t + niJO ] 

m=- 3 n=-“ 

(6) 

where A m , Bm, C m are time-varying, periodic, complex coefficient matrices, and F n 
is a complex vector. These coefficients are shown in appendix C. Let us assume the 
solution of equation (6) as 

00 

Aq = Y q^ exp[i(w G t + ki|i)] (7) 

k=-°0 

where q k is a complex, constant vector. If equation (7) is substituted into equa- 
tion (6) and the same harmonic components on either side are equated, we have an 
infinite number of simultaneous equations involving qk (k = . . . , «) . 

00 00 

Y Y f - ( w G + kfl) 2 ^ + i(w G + k^)B m + C m ]q k exp[i(m + k)ip] exp (iw G t) 

k=-co m=-°° 

00 

' = !C F n exp(in^)exp(iw G t) (8) 

n=-oo 

By eliminating exp(iu> G t) from this equation, equation (8) becomes 

00 00 00 

Y Y Rmk^k exptifa + k)ip] = Y F n ex P( in ^) (9) 

k=— °° m=— co n=-°° 

Equating (m + k) to be n on the left-hand side of equation (9) , and after mathemati- 
cal manipulation, the following expression results: 

n+3 

£ 1,-A - F „ (-<■><“) < 10 > 

k=n- 3 

where 

Rn_k (o> G + k^) 2 A n _ k + i(u> G + kfi)B n _ k + G n _ k (11) 

For the practical calculations, we assume that 

q k = 0 for | k | > K (12) 

Then the following equation will be obtained: 
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(13) 
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• 
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F„ t 
V.NJ 


(14) 


where n is truncated at n = N’ . By solving equation (14), constant amplitudes of 
(2k + 1) generalized coordinates, q^, can be determined, and the blade perturbation 
motion can be approximated by equation (13) . 


Hub Forces and Moments 


In the previous section, the blade equation of motion and its solution using the 
harmonic method are discussed for a sinusoidal vertical gust with frequency o>q. In 
this section, hub forces and moments of a rotor responding to a gust are derived. 
Therefore, we can refer to these responses as "rotor impedance" in the general sense. 
In this study, airloads are obtained by integrating section loading from two- 
dimensional quasi-steady thin airfoil theory. The lift deficiency function C' is 
assumed to be unity. 


The solution of the generalized coordinate Aq is expressed by equations (13) 
and (14) , and perturbations of blade deflections from a trim position are given by 
subtracting an equilibrium value from each expression in equation (1). Let vectors 


AF^ and AM^ denote the variations in hub forces and moments, respectively, owing to 
The augmented vector [AF^ ,Aif£] 


the Jlth blade, 
form in the rotating hub frame: 


can be expressed in the following 



00 

][) {Z 1 }exp[i(w G t + n^)] + {Z 2 )Aq + {Z 3 >Aq + {Z 4 >Aq 

n=-a> 


(15) 


where {Z^ is a complex vector which gives the loads of a completely rigid blade, and 
{Z 2 }, {Z 3 )> and (Z 4 ) denote complex matrices that account for the effects of the 
elastic blade displacement. The terms { Z x } through {Z 4 ) are all functions of the 
azimuthal angle of the Jtth blade, \p £, and each of them has harmonics in U P to 

the third harmonic as follows: 
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3 *\ 
<Zr> = E z im ex P( iim ^> 


m=-3 
3 

= E Z 2 m exp (imp z ) 

m=-3 

3 

{Z 3 } = 2 z 3m exp(innj^) 

m=-3 


{Z 4 > = ]£ Z 4m ex P 

m=-3 


(16) 


J 


Substituting equations (7) and (16) into equation (15) results in the following form: 

{ nr A 3 r °° 

> " E \ E z im exp[i(o) G t + n^ £ + m^)] 

AMj H m=-3 vn=-«> 

+ Jjj [~(w G + kfi ) 2z 2m + 1 (“g + kfi ) z 3m + Z i+m^k ex P[i(w G t + kip z + mij^) 


k=-00 
3 t oo 


3 f oo 

E I E Z im exp{i[w G t + (n + m)^]} 

m=-3 Vn=-°° 


+ E z Sl A exp{i[u> G t 

k=-°o 


+ (k + m) ip « ] } 


} 


(17) 


where 

Z 5m ~^G + k ^) Z 2m + + k ^ Z 3m + Z 4m (1®) 

Here Z im and Z sm denote a constant vector and a constant matrix, respectively. 

The hub loads in the rotating hub frame can be transformed to the rotor frame as 
follows: 



(19) 


where Trr is denoted by equation (A6) in appendix A. Introducing 


gives the following form: 
. 1 3 f o o 

7 - £ £ | £ 

AMjjJ p=-i m=-3 ^n=-« 


KpZ 1Itl exp{i[o) G t + (n + m + p)^]} 




+ ]£) sm^k ex pU[w G t + (k + m + p)^ £ ]} 
k=-» 


i (n+p)+3 “> 

= E ,E, E Vi.v-Cn+p) 6 ^^^ + 

p=-i v=(n+p)-3 n=-°° 


: v = n + m + p 


1 (k+p)+ 3 00 

+ E E E Vs.v-W’k exp[i(w G t + v^)] : v « k + m + p 

p=-i v=(k+p)-3 k=-“ 


( 21 ) 


where Kp is a complete matrix. 

For a practical calculation, the limits of n and k are taken as N' and K, 
respectively. From equation (21), the following equation can be obtained: 


fAF.1 f N'+4 N' 

I V -\ E, E , 

I™*! b*- B n "- H 


K+4 K 

T vn + E E 

v=-K-4 k=-K 




M+4 


E Z v expti(co G t + v^)] 

v=-M-4 


( 22 ) 


where 


i 




T »»- E Vl,v-(n+p) 


P = “l 
1 


s vk E K p Z 5 ,v-(k+p) 


> 


(23) 


p=-i 


M 


( isr * (n' > k) 

K (K > N’) 


where Z v is a vector and denotes the vth harmonic component of the hub loads in 
the rotor frame. Since the Fourier coefficients have been truncated at the N'th or 
Kth harmonic in the blade motion, the hub loads should also be truncated at the 
N'th or Kth harmonic. The total hub loads can be obtained by summing up equation (22) 
for all blades. Using the multiblade summation formula 


8 


N 


£ exp(ivijj £ ) 
SL=i 




exp(ivif)*) : v/N = integer 

0 : otherwise 

the following final expression can be obtained: 


( 24 ) 


N 

E 


= N • 2-r Z N exp[i(w G t + qNij)*) ] 

R k|£M/N qK 


(25) 


where N denotes the number of blades, q is an integer, and if is the azimuthal 
angle of the reference blade. 

Individual Blade Pitch Control and Stability Analysis 

The feedback control system which can control blade pitch individually is dis- 
cussed in this section. Let us consider the following feedback system: 

A0 = (K Aw + K* Aw + K- Aw) + (K Av + K* Av + K~ Av) + (K, Ad) + K* AA + Kv A<j>) 

W W W V V V if> ^ <f> 


= K? Aw(r) + kJ Av(r) + A^>(r) 


(26) 


where 


K = (K , K., K") 
w w w w 


K = (K , K., K-.)' 
V V V V 


K = (K,, K* , K”)- 
9 9 9 9 


Aw(r) = [Aw(r), Aw(r), Aw(r)] 1 
Av(r) = [Av(r) , Av(r), Av(r)] T 
A if) (r) = [A<f>(r), A<£(r) , A$(r)] T 

V Ky, and are feedback gains and Aw(r), Av(r), and A<f>(r) are flap, lead-lag, 
and torsion blade deflections from trim position at r, respectively. By using 
equations (1) and (5), Aw, Av, and A<j) are given as follows: 


Aw = 23 

3 = 1 

(5 j 

Mj. 

W. 

J 

Aq r 

w. 

J 

Aqj) 

00 







av = y] 

j=i 

( b 

Aqj. 

b 

Aq r 

v . 
3 

Aqj) 

S W£ 

h 

-e- 

<i 

( b 

Aqj, 

b 

Aq j , 


Aij) 


( 27 ) 
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For simplicity, let us consider only the flap deflection of a blade; that is, 


Then A0 becomes 


and Aq is 


K = K , = 0 
v <j> 


A0 = K Aw + K Aw + K» Aw 
WWW 


K 


Aq = (Aq x , Aq 2 , . . ., Aq L ) = £ q k exp[i(io G + kft)t] 

k=-K 

From these equations, the following expressions are obtained: 


( 28 ) 


K ^ 

Aw = 53 (w a , + " 2 q ?k + • • • + w L q Lk )exp[i(o> G + kfi)t] 
k=-K 


R 

Aw = i ,(u) G + kfi)(w 1 q lk + . 
k=-K 

K 

Aw = - (m G + k^) 2 (w 1 q lk + 
k=-K 


+ w L q Lk )exp[i(w G + kfi)t] 


• • • 


+ w L q L k> ex p[ ;L ^G + kfi ) t 3 


(29) 


Substituting equation (29) into equation (20) results in the following form: 


A0 = K Aw + K- Aw + K- Aw 

WWW 


K 


= 53 A0 k exp[i(a)Q + k0)t] 

k=-K 


= A0 O exp(iw G t) + A0_ x exp[i(a)Q - fl)t] + A0 X exp[i(w G + ft)t] + . . . (30) 

In equation (30), A0 k (k = -K, . . . , K) are considered as control inputs. There- 
fore, appropriate combinations of A0 k may enable us to alleviate the oscillatory 
gust response of a helicopter rotor. The extension of this idea to an optimal blade 
control algorithm will be discussed in the next section. 

By substituting equations (26) and (27) into equation (C5) (appendix C) , the 
following equation can be obtained: 

53 (A! . Aq . + B’ Aq . + C' Aq.) = F! (31) 

V ij i ±3 ij J ' 
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where 


A'. . 

= A. . • 

- D. 

(K.. 

Aw 

+ 

K- 

Av 

+ 


Ac P) 



J 

w 



V 




b! . 

= B , . ■ 

- D. 

(K. 

Aw 

+ 

K. 

Av 

+ 

K ! 

Ah 

ij 

ij 

1 

w 



V 



<l> 

c! . 

- C. . ■ 

- D. 

(K 

Aw 

+ 

K 

Av 

+ 

k a 

Acj.) 

1 1 

ij 

1 

w 



V 





F' = F. - D. A0 
J J J 


fR 

°j = J {( *j + ab ^j )c aU T + [Wj + (1/2 + a)b<j) j ]c 1 U^ 

- m(k 2 + e Q e)fi 2 ^ + mefi 2 Wj + me6vj}dr; 





+ c x U T U P " 


J 


( 32 ) 


In this feedback system, each sensor (that detects blade deflection) is assumed to be 
located at radial position r along the blade span. Signals from each sensor have 
trim values subtracted in order to get perturbation quantities such as Aw, Av, A<j>. 

By multiplying these perturbation values by feedback gains, this system generates 
inputs to be fed back to the blade pitch control. In this process, any time lags of 
the sensor /actuator are not considered. 

Now let us consider the stability of the blade equation of motion expressed by 
equation (31). In this equation, each coefficient, Ajj , B ij’ and C.-H, Is time-variant 
with the fundamental period T = 2ir. In this study, Floquet theory (refs. ; 97-101) is 
used in order to analyze: the stability of a blade. Define a new variable set as 


X = 



(33) 


By this transformation, equation (31) becomes a first-order differential equation 
with zero external force, as shown by the following: 


[X] = -[E..(t)][X] 


(34) 


where 


[E ij (t)] 


rtAUiB* : 

1J i j 
“I 


[Ai.HOy 




(35) 


and I is a unit matrix. According to the Floquet-L iapunov theory, we can assume 
the following :solution: 

[X(t)] = [G(t)l][a k exp(n k t)] (36) 

The iG(t) is a periodic function with period T, and the n k are complex character- 
istic numbers. The a k are constants derived from initial conditions, that is: 
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(37) 


(a k > - [GCOr^XCO)] 


The Floquet transition matrix [Q] of this system is defined by 

tX(T) ] - [Q] [X(0) ] 


(38) 


The stability of the system, therefore, can be completely determined from the values 
of q k . Inserting the relation for {a k } of equation (37) into equation (36), one 
obtains: 


[X(t)] - [G(t)][- exp (n k t)][G(0)r 1 [X(0).] 

Introducing the state transition matrix [ cj) (t ) ] defined by 

[x(t)] = [<Kt)][x(o)] 

[*(0)] = [i] 

one obtains by comparison of equations (39) and (40) 


(39) 


(40) 


[<f>(t)] = [G(t)][-exp(n k t)][G(0)]- 1 


(41) 


or 


[G(t) ] = [<|)(t)][G(0)][-exp(-n k t)] (42) 

Setting t to be T and considering [G(T)] = [G(0)] gives 

[GWr^CDHGtO)] = [-exp(-n k T)] (43) 

Thus 


[<J> (T) ] = [Q] (t = T) 

Substituting equation (40) into equation (34) gives the following equation: 

d[<j>(t)]/dt + [E ±j (t)][<j>(t)] = [0] 


(44) 

(45) 


Equation (45) is solved by numerically integrating from t = 0 to t = T, with ini- 
tial conditions given by equation (40) . Eigenvalues A k are given, by using equa- 
tions (43) and (44) , as 


[~A k ] = [-exp(-q k T)] - [G(0) ] -1 [Q] [G(0)] 

[Q] = I*(T)] 


(46) 


Numerical Calculations 

In this section, numerical calculations are conducted to determine the capability 
of the individual blade pitch control system. The properties of the blade used in 
this study (an H-34 rotor blade) are given in table 1. The blade structural charac- 
teristics and rotor operating conditions are given in tables 2 and 3, respectively. 
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Gust shape and sensor locations on a blade are schematically depicted in figure 1. 

The hub forces (F x , Fy, F z ) and the hub moments (M x , My, M z ) in the rotor frame are 
shown in figure 2. Blade natural frequencies and mode shapes were calculated using 
the Holzer-Myklestad method (ref. 9). Figure 3 shows the natural frequencies and 
mode shapes for the rotor blade. 

All quantities used in this calculation are expressed as complex numbers, for 
example, Wq = wq g exp(iwQt). The undeformed elastic axis is assumed to lie in the 
X H Z H plane and to coincide with the aerodynamic center axis (e 0 = 0.0, a = -1/2). 
Eight coupled modes (£ = 8) were used (see fig. 2). In figure 3, the flapwise bend- 
ing modes (Nos. 2, 3, 5, and 7) and the lead-lagwise bending modes (Nos. 1,-4, and 8) 
correspond to those of a blade with collective pitch angle of 8°; the sixth mode is 

a torsional mode. The highest harmonic numbers are assumed to be truncated at 5 for 

the gust expression (n = 5) and for the solution of the generalized coordinate 
(k = 5) . Correspondingly, |q| = M/N reduces to q = 0, ±1 in equation (25). Under 
the operating conditions, the equilibrium positions (in this case, the time- 
independent rigid body portion) of the blade deflection are as follows: the flap 

deflection (w^) is 0.11 rad, the lead-lag deflection (Vq) is -0.017 rad, and the tor- 
sional deflection (<J> 0 ) is -0.0004 rad. These values were used in the calculation in 
order to evaluate the effects of the in-plane Coriolis force. The flowchart of the 
computer code for the harmonic method is given in appendix D. 

Figures 4(a) to 4(g) show the rotor response to a sinusoidal gust. The rotor 

impedance is defined as the amplitude ratio between the additional responses of the 
rotor from the steady responses and those owing to sinusoidal gusts with unity ampli- 
tude (wqq = 1.0 m/sec). Phase shift to a vertical sinusoidal gust is also shown in 
figure 4. The horizontal axis denotes nondimensional gust frequencies divided by the 
rotor rotational speed, which varies from 0.0 to 1.0. Each figure (figs. 4(a)-4(g)) 
includes the responses at frequency (u>q) and (wq ±4Q) . The response at frequency 
(wq) is hereinafter called the major response (shown by solid lines in the figures). 
The responses with frequency (wq - 4fi) and (wq + 4Q) are called the subharmonic 
response and the superharmonic response, respectively. In the case of the horizontal 
force F x (fig. 4(a)), each response increases gradually as nondimensional frequency 
(wq/^) increases. The major response has a peak at about 0.7 and has a minimum 
around 0.9. This phenomenon comes from the effect of the lag dampers. Without any 
lag dampers for an articulated blade, the responses will diverge at lag frequency ^ 

(co^ - fi) . The lag dampers suppress this extremely large response. The symbols (V) 
on the horizontal frequency axis denote the location of (w^/fi) and (wl - fi)/fi. The 
subharmonic response is slightly greater than the superharmonic response in the non- 
dimensional frequency range of 0.2 to 0.8. In the case of the side force Fy 
(fig. 4(b)), the tendency of each response is similar to that in the case of the 
horizontal force. There is a sharper reduction in response at (wg/fi) =0.9 for the 
major response. For both figures 4(a) and 4(b), the major response is dominant up 
to 0.8. The subharmonic response is again greater than the superharmonic response 
in the nondimensional frequency range of 0.2 to 0.8. For the vertical force F z 
(fig. 4(c)), a significant difference can be seen relative to the other impedances 
(F x , F ). The major response decreases gradually as the gust frequency increases. 

The subharmonic response is almost the same as the superharmonic response. The 
amplitude of the major response of F z is greater than that of F x or Fy. Quali- 
tatively, the major response of the vertical force is from 2 to 100 times larger than 
the other responses below the nondimensional frequency 0.5. As the frequency 
approaches 1.0, the three response amplitudes become similar. 
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Concerning the phase angle characteristics of the hub forces, the phase shift 
of the subharmonic response is relatively insensitive to the gust excitation com- 
pared with the other responses. From these calculations, the angles are about +90° 
for the horizontal force, about -150° for the side force, and about -90° for the 
vertical force. On the other hand, the phase shift characteristics of the major 
response depends on the gust exciting frequency. Those angles vary from -180° to 
+180° over the frequency range for the horizontal force. However, in the case of 
the vertical force, it seems to be insensitive to the vertical gust. It is almost 
constant at about -10° to 20° below the nondimens ional frequency 0.5. These results 
agree with experimental results (ref. 9). For the case of the side force, those 
shift angles vary from +90° to -120°. 

For the hub moment responses (figs. 4(d) to 4(f)), the trend of the impedance 
is similar to that of the vertical force. Below the nondimensional frequency 0.5, 
the major response dominates the other responses. As with the case of the horizontal 
force and side force, a minimum around the nondimensional frequency 0.9 can be seen 
in the major response of the rolling moment M^.. In the phase angle characteristics 
of the hub moments, the angles depend on the gust exciting frequency. 

In figures 5(a) to 5(c), root loci of the blade characteristic equation with 
feedback (the individual feedback gain is fed back independently) are shown. Floquet 
theory is used in this stability analysis to solve the characteristic equation with 
time-variant, periodic coefficients. The feedback gain is assumed to be a real 
number. Figure 5(a) shows the case in which only the flap deflection (Aw) is fed 
back to the individual blade pitch angle. The feedback gain (i^) varies from 0.0 to 
0.15 rad/m. The limit of the feedback gain is about 0.111, above which the feedback 
system becomes unstable. Figure 5(b) shows the case in which only the velocity of 
the flap deflection (Aw) is fed back. The feedback gain (K^) varies from 0.0 to 
0.005 rad/m/sec. The point at which the system becomes unstable is about 
0.00515 rad/m/sec. In figure 5(c), root loci for the case of the acceleration feed- 
back (Aw) is depicted. The feedback gain (Kft) varies from 0.0 to 5><10 _3 rad/m/sec 2 . 
In this region, the system is stable; for (K^) > 5xl0 -3 the system is unstable. 

Figure 6 shows the sensitivity of the phase angle of the feedback gain to the 
vertical hub force response at gust frequency 0.5 Hz (wq/O = 0.265) and amplitude 
1.8 m/sec. The vertical axis denotes the amplitude of the oscillatory vertical force 
or the nondimensional vertical force (thrust coefficient). The horizontal axis is 
the phase angle of the feedback gain, <p, which varies from 0° to 360°. In this cal- 
culation, the amplitude of the feedback gain is assumed to be constant. The dashed 
line in figure 6 shows the vertical force response without any feedback system (base 
line) . The feedback system used here generates the control input to be fed back to 
the blade pitch according to equation (26). Each sensor is assumed to be located j 
at the blade tip and measures blade deflection independently. Each measured deflec- 
tion (Aw, Av, Acf>) is assumed to be fed back to the blade pitch. The amplitude of 
the vertical force depends on the phase angle, and it has a minimum point at the 
appropriate combination of the amplitude and the phase angle of the feedback gain, 
that is, (K^, (j>) , (Ky, 4>) , etc. From these calculations, the minimum point may be 
realized at <p = 30° for Aw feedback, <f> = 30° for Av feedback, and <j> = 210° 
for A<f) feedback. For A<j> feedback, | j =5.0 was used because the torsional 
deflection of the blade is very small compared with the other deflections. 

Figures 7(a) to 7(c) show the sensitivity of the amplitude of the feedback gain 
to the vibratory vertical force response at the same operating condition. In fig- 
ure 7(a), the phase angle <}> is assumed to be zero; ^ varies from 0.0 to 
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-1.6 rad/m. The amplitude of the vertical force at first rapidly decreases, then 
slowly decreases as further decreases. In this figure, the symbol ® denotes 

the calculation results by the local momentum theory (LMT) . In the case of Av 
feedback (fig. 7(b)), the phase angle (f> is assumed to be 30°, and Ky varies from 
0.0 to -0.1 rad/m. The amplitude of the oscillatory vertical force linearly 
decreases and reaches a minimum at about Ky = -0.04; then, it increases linearly as 
Ky further decreases. For the ■ A<j> feedback (fig. 7(c)), the general trend is simi- 
lar to that in the case of Aw feedback. However, Ka is several orders of magni- 
tude greater from that shown in figure 7(a). This difference is due to the magnitude 
of the blade deflection; that is, usually the torsional deflection is very small com- 
pared with the flap deflection. From these three figures, it is found that the ver- 
tical force response is not as sensitive to the variation of or for Aw or A<|> 

feedback as a variation in Ky for Av feedback. It follows that it is desirable 
to select the Aw or A4> deflection as the measuring signal in order to alleviate the 
vertical hub force response to vertical gusts. 

Figures 8(a) to 8(c) show the effect of the sensor position on the hub response 
to the specified combination of (I^, <j>) , (K v , <j>) , and (K<j, , ((>) . Four sensor positions 
are investigated: r = (1/4)R, r = (1/2)R, r = (3/4)R, and r = R. In these figures, 
the base line is depicted by a dashed line. An effect of sensor position is shown in 
the vertical force and the yawing moment responses. Positions of the sensor farther 
outboard are more effective in the gust alleviation system. It should be noted that 
there is an increase in some responses because of the individual blade pitch control. 
However, this increase is small compared with the vertical force and the yawing 
moment reduction. 

Figures 9 to 14 show the hub gust response calculated by the LMT. Each figure 
includes the three hub forces (F x , F y , F z ), the three moments (M x , My, M z ), and also 
the pitch angle input of a reference blade. The same operating conditions used in 
the harmonic method are used in this calculation.' Gust shapes assumed here are sinu- 
soidal (figs. 9-12) and step (figs. 13 and 14). These gusts are applied to the cal- 

culation after the thrust level has reached a steady-state trim value (usually 
requiring six rotor revolutions for the LMT). The gust amplitude is 1.8 m/sec, which 

is almost half the mean induced velocity of the rotor. The time history of the gust 

shown in each figure is measured at the rotor head. The hub responses are measured 
from when the reference blade is located at zero azimuthal angle. The control pitch 
angle is fed back to the individual blade pitch actuator. Any time lag of the actua- 
tor or other mechanism is not considered. The physical limitation of the pitch 
actuator should be considered for the purpose of implementation. However, in this 
calculation, physical limits were set at the very large values because the emphasis 
of this investigation is on the feasibility of this kind of controller. 

Figure 9 shows the hub gust responses for a vertical sinusoidal’ gust without any 
controller; gust frequency is 0.5 Hz. Among these responses, the vertical force 
(thrust), the rolling moment, and the yawing moment (torque) respond noticeably to 
the vertical gust compared with the others. In the case of the thrust response 
(fig. 9(a)), the amplitude of the fluctuation is about 10% of the steady state 
response. It should be noted that the phase difference between the vertical force 
and the rolling and yawing moments is about 180°. The vertical force responds 
quickly to the vertical gust. 

In figure 10, the hub gust response with the individual blade pitch control is 
shown for the same vertical sinusoidal gust. Feedback gain is -0.5 rad/m. Flap 

deflection (Aw) at the blade tip is selected as a measuring signal. The reduction of 
the oscillatory thrust response is about 60% of that calculated for the uncontrolled 
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case (fig. 9(a)). From figure 7(a), it is expected that the reduction of the thrust 
response should be about 65%. This result implies that the theoretical calculation 
using the harmonic method can reasonably predict the reduction of the thrust response 
In the case of the rolling moment (fig. 10(d)), a significant reduction of the oscil- 
latory response can be obtained. However, from figure 10(f), it should be noted that 
the yawing moment response increases about 20% because of the individual blade pitch 
control. Figure 10(g) shows that the pitch angle of the reference blade varies sinu- 
soidally with a phase difference of 180° relative to the gust. From this result, the 
change of the pitch angle may be expressed as follows: 


0£(t) = [9 0 + A0 Q (t) ] + [0 1C + A0 ic (t)]cos ip + :[e i8 + A0 ls (t)]sin i|T 

A0 o (t) = A0 OC cos(u G t) + A0 Os sin(a)Qt) ( 

A0 lc (t), A0 1S (t) = 0.0 J 

(47) 

where A0 o (t), A0 lc (t), and A0 ls (t) are control inputs generated by the controller, 
and all of them are time-variant. The A0 Oc and A0 OS terms are constant values, 
which are automatically decided once the feedback gain is selected. 


Figures 11(a) to 11(g) show the hub gust response for the vertical sinusoidal 
gust without control. The gust frequency of 2.0 Hz is used. The operating condition 
is the same as in the case shown in figure 9. Compared with figure 9, the amplitude 
of the uncontrolled oscillatory thrust response has decreased noticeably (about 60%). 
This decrease is in agreement with the result calculated by the harmonic method 
(about a 63% decrease; fig. 4(c)). Yet the hub forces and moments still respond 
noticeably to the sinusoidal gust. 

Figures 12(a) to 12(g) show the hub gust responses for this vertical sinusoidal 
gust with the individual blade pitch control. The feedback gain is -0.5 rad/m. In 
this case, there is a noticeable reduction in the thrust, side force, and yawing 
moment responses. However, in the rolling and pitching moment response, increases 
can be seen. Contrary to the case of a 0.5 Hz gust frequency, a slight reduction of 
the yawing moment response can be seen in the case of the 2.0 Hz gust frequency. 


Shown in figures 13(a) to 13(g) are the hub gust responses for a step gust. The 
gust amplitude is 1.8 m/sec. The change in the thrust and in the yawing moment 
responses is significant. The thrust response increases gradually to the new steady 
state without any overshoot. The yawing moment response, however, shows an overshoot 
and then reaches a steady state value. Other responses do not show any significant 
changes. Figures 14(a) to 14(g) show the hub gust responses with the individual 
blade pitch control for this step gust. The feedback gain used in this calcula- 

tion is -0.5 rad/m and the phase angle is assumed to be zero. The general tendencies 
of the response are similar to those of the case shown in figure 13. For the thrust 
response, a reduction of about 70% is achieved. The difference between the steady 
state value and the trim value for yawing moment (fig. 14(f)) is roughly twice that 
shown in figure 13(f). A collective pitch angle change of about 0.6° is input to 
maintain the trimmed condition of the rotor. 
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3. ADAPTIVE BLADE PITCH CONTROL 


In this section, an adaptive blade pitch control designed to alleviate gust- 
induced vibration is analytically investigated for a helicopter with four articulated 
blades. Multicyclic control (MCC) and higher harmonic control (HHC) have previously 
been applied to the reduction of inherent responses of a helicopter, such as N-per- 
revolution vibratory responses. In these controllers, pitch inputs at harmonics of 
the rotor rotational speed are used as control inputs. For the gust response of a 
helicopter, additional vibratory responses, such as gust harmonic (gjq) , subharmonic 
(wq - Q) , and superharmonic (tog + fl), appear in the response. Gust harmonic response 
is the most dominant in the comparatively low gust frequency range (refs. 5 and 9). 

In this study, this dominant response is selected as the response to be reduced. 

An adaptive blade pitch control system is based on digital optimal control 
theory. This system is composed of a set of measurements (oscillatory hub forces), 
a control system based on the optimization of the quadratic performance function, 
and a simulation system of the helicopter rotor. In the following subsection, the 
details will be explained. 

Helicopter Model 

In this study, the helicopter is represented by a linear, quasi-static frequency 
domain model relating the output z to the input 0. Here, z is a vector of the 
harmonic of the gust-induced vibration in the rotor frame. The input 0 is at the 
frequency corresponding to the gust frequency. It is assumed that the gust is sinu- 
soidal, and that the gust frequency is known a priori. When a rotor penetrates into 
the gust, the response at the hub position can be considered to be composed of gust 
and control inputs, such as 

{z} = {z G > + {z c > (48) 

Equation (48) can be rewritten by using the expression of the rotor impedance [Tq] 
and rotor transfer function [T G ] as follows: 

{z G } = [T g ] { g) 

{z c } = [T c ] { 0} 


That is, equation (48) becomes 


where 


{z} = [T G ]{g> + [T c ]{0> 


{z} = (m x l) 

[T g ] - (m x 3) 

{g} = (3 x 1) 

[T c ] = (m x n ) 

{0} = (n x l) 


vector 
matrix 
vector > 
matrix 
vector^ 


(49) 


( 50 ) 
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Once the characteristics of a gust are determined, [Tq ] becomes a constant matrix. 
Therefore, we may denote the first term of equation (49) as {z 0 >. Then equation (49) 
is 

{z} = {z 0 } + [T c ] { 0} (51) 

This form resembles the global model formulation of helicopter vibration (ref. 70). 
The following "local model" can be also taken into account: 

{Az} = [T c ] { A0} (52) 

where Az means z n - z n _ x , and A0 means 0 n - 0 n _ x . 

From theoretical calculations using the harmonic method, the hub response to a 
vertical gust in the rotor frame (Xr, Yr, Zr) can be expressed as follows: 

Z = Z 0 (t) + Z 1 cos(ojQt) + Z 2 sin(ojgt) + Z % cos (tog - 4Q)t + Z^ sin(oo G - 4fi)t 

+ Z 5 cos (o)g + 4fi)t + Z g sin(o)Q + 4ft)t + H.O.T. (53) 

where Z 0 (t) are the inherent responses of a helicopter rotor represented by 

Z Q (t) = Z Q + Z c cos(Nfit) + Zg sin(Nnt) + H.O.T. (54) 

In these responses, gust harmonic responses (Z lt Z 2 ) T are the most dominant in the 
gust response. Oscillatory hub forces and moments at the gust frequency in the rotor 
frame are chosen here as measurements z: 

{z} = (T c , T g , H c , Hg, M yc , Myg, Y c , Y g , M^,, M X g , n zc , M Z g) T (55) 

where 


AT g = Tg cos(a)gt) + Tg sin(u)gt) 
AH g = Hg cos(wgt) + Hg sin(wgt) 
AM yg = M yg cos(u)gt) + M Y g sin(u)gt) 


> 

AYq = Y c cos(wgt) + Yg sin(u)gt) 

AM xg = cos ( w G t ^ + M XS s^C^g^ 
AM zg = M zg cos(ojgt) + M^g s in (wgt ) ^ 


(56) 


Control inputs are selected at the gust frequency u)g, subharmonic frequency (u)g - ft), 
and superharmonic frequency (uig + ft) in the rotor frame as follows: 

A0 = 0 1 , cos (uigt) + 0 2 sin(o)gt) + 0 3 cos(w G - ft)t + 0 4 sin(wg - ft)t 

+ 0 5 cos(w G + ft)t + 0 6 sin(oo G + ft)t (57) 
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Hence the control vehtor 0 has six components: 

{e} = (0 X , e 2 , e 3 , e 4 , e 5 , e 6 ) T (58) 

These control inputs are superimposed on the blade trim pitch input. 


Identification 


Three cases are distinguished for the helicopter model, depending on the identi- 
fication approach: 

1. Identify {z Q } only 

2. Identify [Tq] only 

3. Identify {z Q } and [T^] 

For the gust response, the uncontrolled response {z Q } is generally time-variant ; 

the matrix [Tq] depends on the operating flight conditions. Hence, it is necessary 
that the transfer function [T^] be identified simultaneously with an uncontrolled 
response {z Q }. Case 3 is taken into consideration in this investigation. For the 
local model, it is necessary that the [Tq] matrix be identified for each time-cycle. 
The Kalman filtering technique (refs. 102 and 103) is applied to identify them. For 
the global model, equation (51) can be rearranged as follows: 

z n = z o + ^C®n = [T C , Zo ](0 n ,l) 

= T9 n (59) 


where 


T = m x (n + 1) matrix 


0 n = (n + 1) x 1 vector 

and the subscript n denotes the time-step at t = n At. For this study, it is 
assumed that there is no noise in the measurement of 0 n . The identification algo- 
rithm can be derived by considering the jth measurement as 


z jn 6 n t jn + v jn 


(60) 


where t-j n is the jth row of T, and Vj n is measurement noise, which has zero 
mean, 

E(v n ) = 0 (61) 


and variance 



and a Gaussian probability distribution. Here 6 nm is the Kronecker delta function, 
and the subscript j will be omitted to simplify the notation. The variation of the 
parameters will be modeled as a random process: 


fc n+i fc n 


+ u n 


(63) 


where u n is a random variable, which has zero mean 


and variance 


E(u n ) = 0 


E K%> = Qn 6 


nm 


(64) 


(65) 


and a Gaussian probability distribution. The minimum error- variance estimate of t n 
is then obtained from the Kalman filter (ref. 103), 


where 


£ n = € n-i + k n< z n " 0 n £ n-i> 


«n " Vi + Vi 


P = M - M 0 0 T M /(r + 6 T M 0 ) 

n n nnnn n nnn 


k = M 0 / (r + e„M0 ) 
n n n ' n nnn 


(63) 

(67a) 

(67b) 

(67c) 


Here is the variance of the error in the estimate of t n before the measurement 

and P n is the variance after the measurement; is the Kalman gain vector. To 

simplify the calculation, it will be assumed that Q n and r n have the same time 


variation for all measurements, and that 
same function f j : 

Qn» r n> and p o 

are proportional to the 

r jn 

= f .r 
j n 

(68a) 


- f A 

(68b) 

P jo 

= f i P o 

(68c) 


Then it follows that Pj n = fjP n and Mj n = fjM n ; and that the gain is also 

then the same for all measurements. With the same gains, the rows can be combined 
to form 


T = T + (z - T tO^)^ 
n n-i n n-i n n 

For the local model, equation (69) becomes the following: 


(69) 


T = T + t(z - z ) - T (6 - 6 )]k 

n n-i n n-i n-i n n-i n 


(70) 
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k 
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n n 

^n 

i 

e 

ii 
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Controller 


(71a) 

(71b) 

(71c) 

(71d) 


The control algorithm is based on the minimization of a performance index J 
that is a quadratic function of the input and output variables. The quadratic per- 
formance function to be used here is 


J = z^W z + eJ;W n 0„ + A0^W . q A 0 
nzn n0n nA0n 


(72) 


where W z , Wq, and W^q are weighting matrices, which are assumed to be diagonal and 
the same value for all harmonics of a particular quantity. Then J is a weighted 
sum of the mean squares of the gust response and control. The matrix Wq constrains 
the amplitude of the control, and W^@ constrains the rate of change of the control. 


For the deterministic controller, the control required to alleviate the helicop- 
ter vibration is given by substituting for z n in the performance function J, using 
the helicopter model, and then solving for 0 n that minimizes J. For the global 
model (eq. (51)), the solution can be obtained as follows: 


where 


e n " Cz o + C A0 0 n-i 


C = -DTpW 
o z 


C A0 DW A0 

D - ( T Jw z T c + W 0 + W ^)- 1 

for the local model (eq. (52)), the solution can be obtained as follows: 

0 n = Cz n-i + ( C A0 " CT C ) ®n- l 


(73) 


(74a) 

(74b) 

(74c) 


(75) 


In this derivation, the response z is assumed to be deterministic; therefore, it is 
referred to as the deterministic controller. When the parameter uncertainties are 
taken into account, the cautious controller can be obtained by using the expected 
value of the performance function: 
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J - + 9 X 9 n + 49 n W A6 49 n 

' E j£ <Vi4i] + 9 n w 6 9 n + “5«18 4t n 


where it is assumed that W z is diagonal, and 0 n remains deterministic, 
case of the open- loop control (z Q feedback) , there follows 


(76) 
For the 


.[e s -Iv'i 


L w . (z . + ejt. ) 2 + w . (M + 20-Jk + e^M.. 0 ) (77) 

zj jo n jn zj s zz n tz n tt n- 


where 


M = 


[ 


M M ' 
tt tz 


M . M .. . 
zt zz J 


(78) 


M, 


tt = (n x m) 


matrix 


M tz = M t = (n x 1) vector 


M = scalar 
z z 


So the performance function becomes 


J = z n W z z n + 6 n W 0 e n + ^ (I/ W zj M jnj (j + A6 n W A0 A0 n 
The solution for the control that minimizes J is then 


(79) 


e„ = Cz ~ + c AA 0 „ , + c 

n o A0 n-i o 


(80) 


where the gain matrices C and C^q are the same as for the deterministic controller, 
using the identified values of the parameters and with Wq replaced by 


W fl + £ w . (M. . ) . 

0 zj ' tt y jn 


(81) 


The new constant term is 


c „ - -D [E » 2j 


(82) 
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Similarly, for the case of the closed-loop control (z n _ 1 feedback), the performance 
function is 

J - z *Vn + 9 JV„ + K(»Af> + E (83) 

The solution is identical to that for the deterministic controller, using the identi- 
fied values of the parameters, and with W^q replaced by 

W 46 + E »z 3 M jn <“>> 

j 


Regulators 

A controller combining recursive parameter estimation and linear feedback is 
called a self-tuning regulator. There are two fundamental options for the identifi- 
cation: the use of either an invariable algorithm or an adaptive algorithm. There 

are also two fundamental options for the controller: open-loop or closed-loop. 

Hence, there are four possible regulator configurations (fig. 15). For the adaptive 
algorithm, the parameters are recursively identified on-line, using a Kalman filter. 
For the open-loop algorithm, the control is based on the uncontrolled vibration level 
z Q (identified either on-line or off-line). For the closed-loop algorithm, the con- 
trol is based on the feedback of the measured vibration, z n _x. These four regulator 
options are reviewed in reference 70. 

In this study, two regulators are investigated for a gust alleviation system. 
Figure 16(a) shows an adaptive open-loop regulator algorithm, in which both (z 0 ) and 
[Tq] are identified by a Kalman filter. Figure 16(b) shows an adaptive closed-loop 
regulator algorithm. In this regulator, only [Tq] is updated by a Kalman filter and 
the output vector z n is contaminated by measurement noise, v n . In both regulators, 
the uncertainty of parameters can be taken into consideration by using the variance 
of the error in the estimate of t n before the measurement, M n +x (see eq. (81) or 
(84)). These regulators consist of on-line identification of the parameters and cal- 
culation of the gain matrix, with feedback of the control based on the measured 
responses. It should be noted that the Kalman filter gains can be calculated in the 
time interval between application of 0 n to the helicopter and the measurement of 
the resulting z n . 


Gust Model 

In the past, studies that have dealt with gust suppression systems have used 
gust models, such as the von Kdrmdn model, Dryden model, step model, sinusoidal model, 
and sine-squares model. These models have some advantages as well as some disadvan- 
tages. The von Kdrmdn and the Dryden models are expressions derived from statistical 
techniques in the frequency domain. These expressions are close to the natural tur- 
bulent flows in the sense of the statistics. However, they are not able to show the 
individual flow pattern of turbulence in the time domain. Therefore, these models 
are not suitable for timewise numerical calculation. On the other hand, the step and 
sinusoidal gust models, etc. , are very simple yet different from actual gust shapes. 
These expressions are easy to handle in numerical calculations. By using these 
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gust models, it is possible to obtain the basic characteristics of the plant model 
for the gust response. In this study, the following gust representation is assumed: 


{g} 



= ( u Go’ v Go’ w Go )T exp (iu) G t) 


(85) 


where ug 0 > v Go» an d Wg 0 are horizontal, lateral, and vertical gust amplitudes, 
respectively; and tog is the gust angular velocity. For the step gust, oig = 0. The 
positive direction for the gust velocity is defined in the rotor frame. Most results 
are for a vertical gust only. For this study, the rotor is assumed to have the 
blades at 0°, 90°, 180°, and 270° when the sinusoidal gust is initially encountered 
by the rotor (the gust velocity field is convected past the rotor by the helicopter 
forward speed). As in reference 20, the ratio of the magnitudes of the gust compo- 
nents is assumed to be 

u Go * v Go ’ w Go : 2.0 : 1.8 (86) 

The helicopter is represented mathematically by a simple linear quasi-static fre- 
quency domain model, relating z n and 0 n . Therefore, a simple sinusoidal gust shape 
is selected in order to investigate the feasibility of the gust alleviation 
algorithm. 


Numerical Calculations 

Numerical calculations have been performed to determine the feasibility of the 
adaptive blade control algorithm. This adaptive blade pitch control algorithm con- 
sists of a set of measurements (oscillatory hub forces and moments), an identifica- 
tion system using a Kalman filter, a control system based on the minimization of the 
quadratic performance function, and a simulation system of the helicopter rotor. The 
operating conditions of the rotor are the same as the case of the individual blade 
pitch control (table 3). The initial estimate of the transfer function [Tq] is 
obtained by using the harmonic method program. The detailed derivation of [Tc] is 
given in appendix E. Tables 4 and 5 show the transfer function of the rotor for 
fg = 2.0 Hz and fg = 0.0 Hz (step gust) ? respectively. In the case of the step 
gust, gust frequency is assumed to be zero and the control inputs A0 and the out- 
puts z are assumed as follows: 

A0 = 0 X + ©g cos ip + Q 3 sin ^ 

T (87) 

{z} = (AT 0 , AH 0 , AM Yo , AY 0 , AM Xo , AM Zo r 

where the subscript 0 denotes the mean value. The [Tg] obtained by the harmonic 
method is regarded as the initial estimate in the adaptive blade pitch control. In 
the controller, the [Tg] is updated by the Kalman filter at every time-cycle. In 
this calculation, it is assumed that there is no noise in the measurement of A0. 
However, measurements z are usually contaminated by noise. Measurement noise is 
generated by a random noise generator. Several sets of noise-to-signal ratios were 
studied before these calculations were made. A noise-to-signal ratio of about 0.05 
was used in this study. After trial and error, the initial values of P D , Q 0 in 
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equation (67) were estimated as follows: P Q = [9.8xl0 4 ]N (or N-m), Q q => [9.8xl0 2 ]N 
(or N-m). The variance of the measurement noise, r, is* assumed to be 9.8N (or N-m). 
To keep the noise-to-signal ratio approximately constant, the following relationship 
is taken into account (see ref. 28): 


r n+i r n^ J n+i^ J n^ 


y* V y K y 

min = n = ^max 


( 88 ) 


where J n is the quadratic performance function at the nth time-cycle. The weight- 
ing matrices of the quadratic performance function are taken as 

W z = [l.oxicr 4 ] 

W 0 = [0.0] (89) 

W AQ = [l.OxlO 4 ] 

Since the dimensional values of the measurements z are used, the elements of the 
weighting matrix W z become small. The pitch control 0 n are unconstrained; the 
pitch control rates A0 n are constrained. The values given by equation , (89 ) were 
selected for good algorithm convergence characteristics. 

In this gust alleviation system, the gust frequency f G is specified. There- 
fore, each time-cycle that updates the parameters must vary depending on the gust 
frequency. In most calculations, a 2.0 Hz gust frequency is used. This frequency is 
more than half of the rotor rotational speed, and somewhat larger than the typical 
mean atmospheric turbulence (usually below 1.0 Hz). An updated parameter estimate is 
performed every four rotor revolutions. During that time, the measurements z are 
discretely sampled by the measuring system. Since the measurements z are data in 
the time domain, it is required that they be converted from the time domain to the 
frequency domain using the fast Fourier transform (FFT) (ref. 104). In the FFT, the 
solution of the data in the frequency domain depends on the sampling time At 
(Shannon's theorem); in this calculation, fiAt = 10°. Therefore, when the gust fre- 
quency f G falls between two points (that is, nAf < f G < (n 4- l)Af, n = integer, 

Af = frequency step) , it is impossible to calculate the correct values by the FFT. 

In this analysis, correct values are approximated from known gust frequency f G and 
by using a cubic spline interpolation method. Furthermore, it is necessary to con- 
sider the calculation time of the optimal control, which determines the optimal con- 
trol for the next time-cycle using the output data from the FFT. One time-cycle is 
8it rad (fig. 17). The sampling interval is (128/144) (8tt) rad. Data are sampled at 
every 10° (144 samples in four revolutions); 128 samples are used in FFT since powers 
of 2 are typically computationally more efficient (2 7 = 128). Consequently the 
required calculation time is assumed to be (16 /144) (8ir) rad. As a result, there is 
a phase shift (<(>j + cpxi ) in the measured response at the beginning of the next time- 
cycle relative to the phase at the end of the previous sampling interval. In the 
program code, this phase shift is accounted for in implementing the control as 
follows: 
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( z i)n = ( z i)n-i cos [ (ojq / f2) 4>n 1 + ( z 2 ) n -i sin[ (w g /Q)(|) II ] 
( z 2 )n = ( Z 2 )n-i cos [ (wg/^)4>h] " (Z^n-i sin[ (u) G /ft) <t>u ] 
Z n-i = ( z 0 ) n -i + ( z i)n-i cos(u G t) + (Z 2 ) n _ 1 sin(u) G t) + 
Z n = (Z 0 )n + (Z x )n cos(m G t) + (Z 2 ) n sin(w G t) + . . . 




( 90 ) 


where Z x , Z 2 are coefficients of the cosine and sine term of the response, respec- 
tively. The subscripts n denote the nth time-cycle. If this phase shift is not 
considered, (Z 1 ) n and (Z 2 ) n coincide with (Z^^ and (Z 2 ) n _ 1 . By using the above 
equations, the optimal control inputs for the next time-cycle can be determined. A 
global helicopter model is used. For the open- loop controller, both z 0 and [T c ] are 
updated by the Kalman filter. Only [T c ] is updated for the closed-loop controller. 


In figures 18(a) to 18(g) and figures 19(a) to 19(g), the hub gust responses 
without control are shown separately for the horizontal and the lateral components 
of the gust. The amplitude of the horizontal gust is 3.0 m/sec, and that of the lat- 
eral gust is 2.0 m/sec. The gust frequency of 2.0 Hz is used. Figures 18 and 19 
show little influence on the thrust response. In the horizontal and side forces and 
yawing moment responses, the effect of these types of gusts on the hub response can 
be seen. Figures 20(a) to 20(g) show the hub responses without control for the three 
dimensional gust. The amplitude of the three gust components are u Go =3.0 m/sec; 
v Go = 2.0 m/sec; and wco = 1.8 m/sec. The three components of the gust have the 
same frequency (f G = 2.0 Hz). The increase of the horizontal force, the side force, 
and the yawing moment responses can be seen by comparing figure 20 with figure 11. 
From these results, a vertical component of the gust is the most influential on the 
rotor gust response (specifically in the thrust response). 


Figures 21(a) to 21(g) show the hub gust responses with the adaptive blade pitch 
control for a vertical gust component only (f G = 2.0 Hz). The controller is deter- 
ministic, and the global helicopter model was used. Two inputs (0 X and 0 2 ) and two 
outputs (cosine and sine elements of the thrust response) are considered. Therefore, 
the dimension of [T G ] becomes (2x2). The thrust response (fig. 21(a)) gradually 
decreases, but the yawing moment response (fig. 21(f)) responds dramatically (see 
fig. 11(f)). The other responses (F x , F y ^ M^, My) are similar to that shown in 
figure 11. In this case, the aim of the controller is to reduce the thrust fluctua- 
tion owing to the gust, and no attempt is made to reduce the other responses. 


In figures 22(a) to 22(g), the hub gust responses are shown with adaptive blade 
pitch control, again for the vertical gust component only. In this case, the dimen- 
sion of [T c ] is (6 x 12). As control inputs, six elements (0 X to 0 g ) are used to 
alleviate the responses. Measurements z involve the hub forces (F x , F y , F z ) and 
the hub moments (M x , M y , M z ). Compared with figure 21, the reduction of the thrust 
response is significant (about 50% to 80%). From this figure, it is observed that 
the fluctuation of the thrust does not reduce uniformly; in other words, it sometimes 
diverges and converges. This phenomenon is explained as follows. The control inputs 
generated by the controller depend on the measurements. When the hub responses are 
decreased by means of the controller, the optimal control inputs necessarily become 
small. However, the gust itself is unchangeable. Therefore, the effect of the con- 
troller on the reduction of the responses becomes weak. At the same time, the uncer- 
tainty of the identified parameters may increase, even if the noise to signal ratio 
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is kept constant by equation (88). In the yawing moment response (fig. 22(f)), the 
characteristics of the response are improved (see fig. 21(f)), and the amplitude of 
the fluctuation is a little less than that of the case shown in figure 11. Similar 
to the thrust response, the response sometimes diverges and converges. This is due 
to the change of the induced drag which is directly related to the thrust. Contrary 
to expectations, the responses (F x , F y , My) show a slight increase in magnitude. 

In the pitch angle change of the reference blade (fig. 22(g)), the sinusoidal change 
with the gust frequency can be observed, similar to that shown in figure 21(g). 

Figure 23 shows the time histories of the mean square thrust response. In this 
figure, the solid line denotes uncontrolled response. The dashed line corresponds 
to figure 21 (two inputs and two outputs) and the broken line corresponds to fig- 
ure 22 (six inputs and 12 outputs). As explained before, the parameters are updated 
at every four revolutions. During one time-cycle, the control inputs are kept con- 
stant. The response for the case of thrust control only is at first constant, and 
then gradually decreases. The response is reduced by almost 50% after 20 revolutions 
with the controller engaged (five time-cycles). The response for the case of complete 
hub response control decreases to a much lower value, although not continuously. A 
50% to 80% reduction of the thrust response is attained for the case of complete hub 
response control. From these results, it is found that the measurements z should 
include not only the thrust response but also other hub responses. 

Shown in figures 24(a) to 24(g) are the hub responses for the vertical step gust. 
Operating conditions are the same as in the case shown in figure 13. The adaptive 
closed-loop controller is used. The measurements z involve only the longitudinal 
response of the rotor; that is, the horizontal force F x , the thrust F z , and the 
pitching moment M y . The control inputs have three elements (0 X , 0 2 , 0 3 ) according 
to equation (87). All parameters (z Q and [T^]) are updated at every rotor revolution 
by the Kalman filter. Significant reduction in 1;he thrust response (fig. 24(a)) is 
achieved by the controller (almost 100%). However, compared with figure 14, the 
horizontal and side force responses transfer to a new steady state condition. The 
rolling and pitching moment responses (figs. 24(d) and 24(e)) transfer to new steady 
state conditions. In the yawing moment response (fig. 24(f)), there is a significant 
change in the steady state condition. This phenomenon is a result of the strong 
effect of the controller on the gust alleviation system. The controller generates 
the optimal control inputs in order to reduce the change from the steady value. If 
these control Inputs are large values, a new trim condition could result. Referring 
to figure 24(g), the maximum change of the pitch angle is about 2°. The cyclic pitch 
angles for the trim condition itself are of the order of several degrees. Therefore, 
the pitch angle change by the controller has an effect on the rotor trim condition. 

In this case, the ratio between the control angles A0 and the trim pitch angle 0 
becomes more than 10%. The flight trim condition consequently changes and transfers 
to the new steady state condition. 

In figures 25(a) to 25(g), the hub responses for the vertical step gust are 
shown. The adaptive closed-loop controller is used, and the operating conditions 
are the same as those shown in figure 24. To investigate the sensitivity of the 
transfer function to the controller performance, arbitrary small initial values of 
the transfer function are used in this calculation. Compared with figure 24, the 
thrust fluctuation (fig. 25(a)) decreases very slowly. For the other responses 
except yawing moments the trim condition remains the same. In the case of the yaw- 
ing moment (fig. 25(f)), the transfer to the steady state can be seen after showing 
an overshoot. This phenomenon is due to the induced drag, as mentioned before. 

From these results, it can be concluded that the performance of the controller 
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depends on the initial estimate of the rotor transfer function. If a more correct 
transfer function is used, convergence of the responses can be obtained more quickly. 
As shown previously, the transfer function derived by the harmonic method gives a 
reasonable initial estimate. 

Figures 26(a) to 26(g) show the hub responses, again for the vertical step gust. 
In this case, the adaptive open-loop controller is used. The measurements z involve 
six responses of the rotor; that is, three hub forces (F x , F y , F z ) and three hub 
moments (M x , M y , M z ). Similar to the case shown in figure 24, all elements of [Tq ] 
are updated at every rotor revolution. The magnitude of the controller inputs is 
constrained by a prescribed maximum value (A0 max =2.0°). Three components (0 lf 0 2 , 
0 3 ) of the control inputs (eq. (87)) are used in this calculation. In this thrust 
response (fig. 26(a)), the deviation from the steady value of the thrust gradually 
decreases after responding to the step gust. Compared with the case of figure 24, 
the pitch angle change is very slow and small. The other hub forces, together with 
the rolling and pitching moments, remain steady. The yawing response deviates slowly 
from the steady value because of the thrust response. In the case of the open-loop 
controller, the optimal control inputs depend on the uncontrolled response z 0 . If 
the z Q are small, then the optimal control inputs may be small. For this case, the 
optimal control inputs at the beginning stage of controller operation are less than 
0.5°. 

In figure 27, the time histories of the quadratic performance function J are 
shown. The solid line denotes the function J without control. The J's with con- 
trol are shown for both the adaptive closed-loop control (corresponding to fig. 24) 
and the adaptive open-loop control (corresponding to fig. 26). For the adaptive 
closed-loop controller, the J decreases rapidly after showing a sharp increase. 

The J for the adaptive open-loop controller decreases moderately after showing the 
same increase. Both control schemes show the effect the controller has on the reduc- 
tion of the hub gust responses. 

Since deterministic controllers have been considered in this analysis , the 
uncertainty of the parameter identification has not been considered in the calcula- 
tion of the optimal control inputs. These cautious properties can play an important 
role for inherent helicopter vibration reduction schemes (ref. 68). In figures 28(a) 
to 28(g), the hub gust responses with the cautious controller are shown for the ver- 
tical gust (fg = 2.0 Hz). Referring to figure 22, there is no significant difference 
between these figures. The cautious controller involves the variance of the error 
in the estimate of the transfer function f T q ] before the measurement. The more 
uncertain the parameter estimates, the larger the variance of the error. In this 
calculation, the term concerned with the cautious properties 

? w zj^tt^jxij 

m 

in equation (81) was about 1% of the term TqW z T(-, in equation (74c). Therefore, 
these calculated results show that the cautious properties have little effect on the 
controller performance in this study. 

Results from applying adaptive blade pitch control for the three dimensional 
gust are shown in figures 29(a) to 29(g). Operating conditions are the same as those 
in the case of figure 20. The large reduction in the thrust response can be observed 
by comparing figure 28 with figure 20. The yawing moment response shows the same 
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characteristics as for a vertical gust only (fig. 22). In general, all the hub 
responses are similar to those of the case shown in figure 22. In this calculation, 
the cautious properties have not been considered because of their minimal effect on 
controller performance. 


4. CONCLUSIONS 


Control systems to reduce the gust-induced vibration of helicopter rotors have 
been studied by two means. One is an individual blade pitch control scheme, and the 
other is an adaptive blade pitch control scheme. Among three gust components, the 
horizontal and lateral components of gust have little influence on the hub gust 
response. However, the vertical component of the gust has a great influence on the 
rotor hub response. In studying the gust response of the rotor, the thrust response 
shows the most significant change compared with the other responses. Hence, reducing 
the thrust response should be the aim of the gust alleviation system. Both control 
schemes show significant effects of the controllers on the reduction of the gust- 
induced vibration of the helicopter rotor. From these theoretical analyses, the 
following conclusions are drawn. 

For the individual blade pitch control scheme: 

1. A computer code that calculates the gust response of a rotor has been 
developed by using the harmonic method. This code is applicable to the rotor impe- 
dance calculation for gust response or the transfer function calculation of the rotor. 

2. Outboard measurement sensors are more effective in reducing the gust-induced 
vibration than are inboard radial station locatipns. 

3. The gust-induced thrust response determined by the harmonic method is about 
10% to 20% larger than that given by the LMI. However, the results by the harmonic 
method can qualitatively predict the characteristics of the gust response. 

4. Individual blade pitch control can significantly reduce gust-induced vibra- 
tion. For example, reductions of the thrust response for the sinusoidal gust by 
almost 60% and for the step gust by almost 70% can be obtained for a specified feed- 
back gain 1^. 

5. For the gust models used, an increase in the yawing moment fluctuation is 
observed. A control scheme to include the reduction of the transient yawing response 
should be incorporated. 

6. The magnitude of the control inputs is less than 1.0°. For these pitch 
angles, there is little effect on the trim Condition of the rotor. 

7. Because the flapwise blade deflection is most sensitive to the gust (com- 
pared with the lead-lag and torsion deflections), it is a good measurement of the 
gust-induced rotor response for an articulated rotor. 

For the adaptive blade pitch control scheme: 

1. The major, the subharmonic, and the superharmonic inputs are considered in 
the controller. As the gust frequency increases, higher frequency terms of the gust 
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response become large. To alleviate these high frequency terms in the vibratory 
response , the number of terms included in the control inputs increases , making the 
controller more complex. 

2. The performance of the controller depends on the initial estimate of the 
rotor transfer function. When the exact transfer function is used, the convergence 
of the hub responses can be achieved quickly. The transfer function derived by the 
harmonic method gives a good estimate. 

3. In using the FFT to convert measurements from the time domain to the fre- 
quency domain, some approximations must be made because of the arbitrary gust fre- 
quency. This may increase the uncertainty of the measurement. A more accurate 
frequency domain determination method is required. 

4. For the case of a sinusoidal gust, the adaptive open-loop regulator is best 
suited for the gust alleviation system. Results show that a 50% to 80% reduction in 
the thrust response can be obtained. The regulator studied in this report is shown 
to be applicable to a three dimensional gust. 

5. For the step gust, the adaptive closed-loop regulator performs better than 
the adaptive open-loop regulator. The closed-loop regulator yields a rapid reduction 
of the gust-induced thrust response (almost 100%), even though it violates the trim 
condition. The open-loop regulator shows that convergence of the thrust response is 
slow. 

6. The uncertainty of the parameter identification has little influence on the 
improvement of the regulator in this investigation. 

7. In this adaptive control scheme, the. gust frequency is prescribed at the 
initial stage of the calculations. Therefore, this type of regulator does not apply 
to random gust responses. 
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APPENDIX A 


COORDINATE SYSTEMS 


Five Cartesian frames are used in this study to describe the blade and hub 
motions. The definition of these frames is briefly reviewed in this appendix. 

Inertial frame (Xj, Yp, Zp)— This frame is inertially fixed and coincides with 
the body frame under the steady flight condition. 


Body frame (Xg, Yg, Zg)— This frame is fixed to the fuselage with its origin at 
the center of gravity; Xg is oriented to the 'direction of the steady flight, Yg is 
taken to be positive in the right-hand direction, Zg is positive downward. The 
relationship between the inertial frame and the body frame is as shown: 
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(A2) 


where (4>, 0 , ¥) are Euler angles of the body, and (Xg 0 , Yg 0 , Zg Q ) is the position of 
the origin of the body frame in the inertial frame. 


Rotor frame (Xr, Yr, Zr)— This frame is fixed to the fuselage with its origin at 
hub center; Zr coincides with the rotor shaft and Xr is positive rearward. This 
frame is related to the body frame by the following transformation: 



-cos;ig, 0, -sin ig 





0 , 1 , 0 

sin ig, 0, -cos ig 


(A3) 


(A4) 


where ig is the shaft inclination angle and(X^ Q , Y^ Q , Z^ 0 ) is the position of the 
origin of the rotor frame in the body frame. 
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Rotating hub frame (X H , Y H , Zjj)— This frame is rotated about the Zr axis by 
ip = fit! The coordinate transformation between the rotor frame and the hub frame is 
expressed as 



(A5) 


cos -sin ip, 0 
sin \p, cos ip, 0 
0 , 0,1 


(A6) 


Blade frame (X]-,, Y-j,, Zfo)— This is a local frame fixed to each blade section con- 
sidered! The blade section itself is assumed to be rigid; Z^ is directed outward 
along the local elastic axis, and Y^, is directed toward the zero-lift angle of 
the blade section. The displacements of the elastic axis are denoted by (u, v, w) 
in the hub frame. Angular changes related to the hub frame are given by the Euler 
angles (0 + <p, -w', v') where ( )' denotes d( )/dXR. The transformation between the 
hub frame and the blade frame is given as 
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+ u"> 

Y " 

► - HbH)" 

V 

► + 

e o + v ► 

LJ 


U> 


l W ; 


[T 


bH J 


1 , -v', 0 1,0, -w' 

v\ 1,0 0,1, 0 

0 , 0,1 w' , 0 , 1 

1 , -v' - 0w' , -w ? + 0v' 

v* , i - o 2 /2 - e<j>, -(e + <j>) 

w' , (0 + 4>) , l - e 2 /2 - e<|> 


1 , 0 , 0 

0, cos(0 + 4>) , -sin(0 + (J>) 

0, sin(0 + <p) , cos(9 + <(>) 


(A7) 


(A8) 


In figures 30(a) and 30(b), four coordinate systems (except the inertial frame) are 
shown schematically. 


Since motion of a fuselage is not considered ($ = 
tion angle is zero (ig = 0), three coordinate systems, 
and (X],, Yfc, Z^) are used in this study. 


0 = Y = 0) , and shaft inclina- 

(X R , y r , z r ), (x h , y h , z h ), 
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APPENDIX B 
BLADE SECTION LOADS 


In this appendix, total loads on the blade section are defined (fig. 31). 

Inertial loads— The inertial forces {P} and moments {Q} acting on a unit span of 
the blade are given as follows: 




{P} = < 


\PzJ 


(Bl) 


'ii - e(v’ + 0w') + 2fi[-v + e0(9 + <j>)] - fi 2 [r + u - e(v' + 0w')] 

Py ^ = -m^.v - e0(0 + $) + 2fi[u - e(v' + 0w’)] - £2 2 [e Q + e + v - e0(j> ] (B2) 

.w + e(0 + $) 


(Q> = 



(B3) 






0 + (|>‘ + fi 2 (0 + (())+ 2fi0v r 
y = -mk 2 ^ -0v r + n 2 0v' 

v' + 0w' - 2Q0 (0 + 40 
/'w - 0[v + 2£2u - fi 2 (e Q + v)] + e^O 2 ^ 

0 [ii - 2 fiv - fi 2 (r + u) ] 

^-[ii - 2flv - fi 2 (r + u) ] - rQ 2 0tf> + e 0 fi 2 (v' + 0w')jj 


me< 


rQ 2 4> 


where 


m = 


l 


dm , 


chord 


if 


n dm 


chord 




(B4) 


(B5) 


are mass, distance between mass center and elastic axis, and polar radius of gyration 
of the blade per unit span, respectively. Equations (B2) through (B4) are the iner- 
tial loads owing .to the blade motion when the hub is not in motion. 

Aerodynamic loads— Aerodynamic forces and moments are derived from two- 
dimensional quasi-steady thin airfoil theory (ref. 96) and assuming airfoil plunging 
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to be zero (h = 0). Blade lift, pitching moment, and profile drag are given as 
follows : 

L = 2TrpbU a C’[H a + aU a + (1/2 - a)b(a + Q 0 ) ] + irpb 2 (d/dt) [H a + aU a - ab(a + Q q )] 

(B6) 

M = (a + l/2)2Trpb 2 U a C' [H a + aU a + (1/2 - a)b(ci + Q 0 ) ] - (Trpb 3 /2)H a (a + Q 0 ) 

+ irpb 3 (d/dt)[a(H a + aU a ) - b (1/8 + a 2 )(ct + Q)] (B7) 

D - pbu|C do (B8) 


D ± = -L sin <j> = ~LH a /U a (B9) 

where C' is the lift deficiency function and is assumed to be unity in this study. 
The quantities D^, <j>, and C^o are induced drag, inflow angle, and drag coefficient, 
respectively. It is noted that in equations (B6) through (B8) differentiation (d/dt) 
must be operated on U a as well as H a , Q 0 , and a. Quantities of U a , H a , Q 0 , and a 
are expressed as follows: 


U a = r£2 + v + U 0 sin ip + v'U 0 cos ip + uq sin \p + Vq cos \p^ 

H a = -Vg - w + wq - w'U 0 cos ip 

■ 

Q 0 = fiw' 

ot = 0 + <j> 


(BIO) 


By using equations (B6) through (B9) , aerodynamic forces and moments acting on the 
blade are 


^ax^ 


■ay 


> = [A A ] 




ax 


1 , 0 , 0 | 

0, cos (p, sin <)> 

0, -sin <p, cos <f> 


0 

l-Dl 

L 


ay 


x aZ' 




r m > 

1! 

O 

► = < 

v'M 

J bJ 

. 

U'M, 


-w'L 
-D - D, 


(Bll) 


(B12) 



1 , -v ' , 0 


> 

o 

t— ! 


> 

*» 

*> 

1 

1 — I 

[A a ] = 

V, 1,0 


0 , 1, 0 


v', 1 , 0 


0,0,1 


w' , 0 , 1 


w* , 0 , 1 


(B13) 
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The wing motions in the two-dimensional thin airfoil theory and the lift, 
pitching moment acting on the wing section are shown in figures 32 and 33 
respectively. 

Total loads— Total loads acting on the blade section can be obtained 
of inertial, aerodynamic, and gravitational loads as 



where g is the gravitational acceleration. 


drag, and 
as the sum 


(B14) 


(B15) 
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APPENDIX C 


COEFFICIENTS OF EQUATION (6): A^, B m , C m , AND F n 

The equation given from equation (2) is 

X) (A 1:j Aq ± + B ± j Aq ± + C^ Aq ± ) = Fj 


(Cl) 


where 


(C2) 


/• R _ . k 2 _ \ /• R 

A ij = m j 5 ij + J VT" C 2^i) dr ~ J (Wj + ab<J>pc 2 (0v ± - v ± - ab$ i )dr 

r R 

B. . = I 2mft [ . (k 2 0v' - e0u.) + w'w.v. + v.{u. + e(v! + 0w!)}ldr 
ij J Q J i i o j i j x x x' JJ 

+ / \ C 2 {Vi + J ^i} dr - jf jwj + (a + |) b^C^OU^, + U p )v ± - U T w. 

+ (t ~ a ) bU T^i + (y ~ a ) b (® + fiw o>^i}dr “ f (w. + ab<j> .)C 2 {U R (0v| - wp 

/ R fR / 

v.C-U^v. dr + I v.C-.uOtL, + 2U )w. - 0U v. 
j 3 Tx J j 1 \ T p i pi 


- (-| - a)bU p 4> 1 + (y - a)b(0 + ^w^w^dr + C ? v|(r ? )vj(r ? ) (C3) 

R *R 

+ 1 (f c ^vi) dr - p{”j + (i + a KK[{- u T u R + (I - a ) baD i} 5 i 

+ (20U t + U )U R vj! + (| - a)b(0 + flw’)U R vj_ + Up^ldr - f (w + ab£ )C 2 

P J * Q J J 

-R 

X {9U R vJ - U e 5’ + 9U r 7! + U T i.}dr + J v.C^-CeU,, + 2D p )li R 5' + (f - a)bU p SK’ 


c. . 

1J 


+ U T U p ^. + 0U p U R v 


*R 

i - (i - a ) b( ® + + J V c 3 u T Vi )dr 


(C4) 
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F. 

J 




ab^)C 2 [-w G ]dr 4- ^ + '{j ~ aJb^C^ 


+ v c,(eu + 2u ) 
3 1 T p 


i (y " a)b(9 + [-w G ]dr + |j [(Wj + ab4>j )C 2 U T + {w.. + + aJbfj^uJ 

+ CjU^U^v., - m(k 2 + e^e)fi 2 |. + mefi 2 w.. + me'Ov^. ] drj A0 + C r v! (r r )v' (r r ) (C5) 


^ T p j 


o 3 


c r v x v 


and where 


C x = 2iTpbC' 
C 2 = irpb 2 


C 3 “ 2pbC do 

b = C/2 

U„ = U cos ip 
R o r 

U,p = rfl + U q sin ip 

U = W - v 
poo 

0 = 0 Q + 8 cos ip + 0 1S sin ip 


Here wq is a vertical gust and A0 is the additional control input to be fed back 
to the actuator. These equations can be expanded, using Euler's formula, to yield 

cos ^ = exp (iiji) + exp (-i;. ) > sin ^ 

where 

i = /=r 


When the two-dimensional quasi-steady thin airfoil theory is used,. there appear har- 
monics of ip up to the third in each coefficient; that is, 

A.. = A exp(-i3ip) + A_ exp(-i2i|>) + A_ exp(-i^) + A + A exp(ii^) 
lj —3 ” 2 *"1 0 1 

+ A 2 exp(i2ip) + A 3 exp(i34») 

3 

= 2 A m exp(im^) (C6) 

m=-3 
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B -H /V 22 exp(im^) 
J m=-3 


(C7) 


c ± a - 22 C m e»P (i>#) 

J m=- 3 


(C8) 


Each term of these coefficients has the ith and jth mode of the blade deflection; 
for example: 


r R - r K - _ 

\_ 1 = -2 J mea s Vj^ dr • 6^ - J (w^ + ab(j)j )C 2 a g v 1 dr 


(C9) 


A o 


r 


m{(k 2 *. + ew. - ea.v.H. + (w, + e<L)w, + (v, - ea e 5.)v.}dr • 6.. 
j j 6 j j j j j j 6 J J ij 

U 

+ / (1 b2(] 2 ^i)^j dr " / ( a g^i " w 1 - ab^ ± )C 2 (wj + ab$j)dr 


(CIO) 


r R - r R 

= -2 J mea 7 v^<j)j dr • - J (w.j + ab<J>j dr 


(Cll) 


A =A =A =A =0 
-3-223 


(C12) 


where 

a 5 = ~2 (®ic + i0 is) 

a 6 = 0„ + 6 t (x - 0.75) 

“7 ' J < 9 1C. * 19 1 S > 

r 

X = R 

and 6,0 t are the Kronecker delta function and pre-twist angle of a blade, respec- 
tively. For the practical calculation, both i and j are truncated by L. Then the 
following equation can be obtained: 


Y) {[A exp (imijO ] Aq + [B exp (iimjj) ] Aq + [C exp(im^)]} = Yj F_ exp[i(u>^t + n^)] 
m=-3 


n=-oo 


(C13) 


where B m , C m are (L x L) matrices in which all elements are complex. 
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Now, consider the expression of a sinusoidal gust in the rotating hub frame at 
the position of hub center. In the rotor frame, a sinusoidal gust is assumed to be 
expressed as follows: 

w G “ w Go exp (iw G t) (C14) 

Then, it can be transformed at position (r,i^) in the rotating hub frame as 

w G = w Gq exp(iw G t) * f(x,^) (C15) 

where 

f(x,ijO = exp[-i(w G /fi) (x/y)cos t|>] 

00 

= E (-i) n J • x/y)exp(in^) (C16) 

n==-oo 

In this derivation, the following formula is used: 

00 

exp(iz sin 0) = E J (z)exp(in0) 

n=-oo 

where J n is nth order Bessel function. By substituting equation (C16) into (C15), 
the expression of the coefficient F n is given by 

V = w Go ( “ 1)nJ n ( V fi -' x/y) (C17) 
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APPENDIX D 


FLOWCHART OF THE HARMONIC METHOD 


The flowchart of the computer code by the harmonic method is shown in figure 34. 
In this program, the feedback gains are specified in advance and read as data. It is 
possible to analyze the stability of a blade and also to calculate the rotor trans- 
fer function using this code. Approximately 1.0 sec of CPU time was required for 
one execution of this code on the CDC 7600 computer. 
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APPENDIX E 


TRANSFER FUNCTION OF A ROTOR 


In this appendix, the derivation of the transfer function of a rotor is 
described by using the harmonic method. First, gust amplitude is assumed to be zero 
in the computer code. Second, the transfer function is assumed to be expressed in 
the frequency domain as follows: 

{z} = [T c ]{e} (El) 

where {z} and {0} are the response vector and the input vector, respectively; {z} and 
{0} are given by equations (55) and (58), respectively. Therefore, [Tg] becomes a 
(12 x 6) matrix. 

In order to calculate the transfer function of a rotor, it is necessary to modify 
equation (31). Let us consider the following control input A0: 

A0 = (Gj - i0 2 )exp(ia)gt) + (0 g - i0 4 )exp[i(u> G - fi)t] + (6 g - i0 g )exp[i(to G + ft)t] 

= { (6 x cos ojgt + 0 2 sin Wgt) + i (0 x sin a> G t - 0 2 cos cOgt)} 

+ {[0 3 cos(a> G - fl)t + 0^ sin(a) G - fi)t] + i[0 3 sin(a> G - fi)t - 0 4 cos(u) G - ft)t]} 

+ [0 g cos (nig + ft)t + 0 g sin(w G + fi)t] + ± [ 0 5 sin(w G + fi)t - 0 g cos(to G + fi)t]} 

(E2) 

The real part of the input A0 corresponds with equation (57). On the other hand, 
the response vector {z} can be obtained by equation (25). From assumptions used here 
(i.e., K = N’ =5), q becomes zero in equation (25). For simplicity, let us con- 
sider only the vertical response: 

AT = Z T • exp (iajgt ) 

= [Z TC cos(o) G t) - Z Tg sin(u> G t)] + i[Z TC sinfagt) + Z Tg cosfagt)] (E3) 

where is the complex response vector. The real part of the is Z^q and 

the imaginary part is Zjg. The comparison of the real part of both response and 
input gives the following relationship: 


(E4) 


where t-y (i = 1, . . ., 12, j = 1, . . .,6) are elements of the [Tq] matrix. Simi- 
lar relationships can be given for the other responses. 


Z TC = 1 1 6 1 + t 12 0 2 + t 13 9 3 + t 14 0 4’ + t 15 0 5 + * 1 *% 
_Z TS = t 21 9 l + t 2 2 6 2 + t 2 3 9 3 + t 24 9 4 + t 25 6 5 + t 26 S 6 
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To determine all the elements of [T C 1 matrix, the least-squares method is used. 
Let us consider the sum of the squares of errors: 


S = 


N 

E 

n—l 


(z . 

jn 


6 nV 2 ■ 


e t.) T (z. - 0 1 .) 
J 3 J 


(E5) 


where N denotes the number of parameters (the dimension of Zj). The number of 
parameters to be identified (the dimension of tj) is assumed to be L. N must be 
greater than or equal to L. Here, the vector Zj and matrix 0 are defined as 



The solution that minimizes S is the least-squares estimate: 

t. = ( 0 T 0 ) -1 G T z. (E7a) 

J J 

or 

tj = z* 0 (0 T 0 ) _1 (E7b) 

Putting the rows together again gives 

T = Z 0 (0 T 0 ) _1 (E8) 


where 



(E9) 


In this study, both N and L have values of 12. Therefore, equation (E8) can 
be rewritten as follows: 

T = Z( 0 T ) _1 (E10) 

Elements for the transfer function [Tq], calculated using equation (E10), are shown 
in tables 4 and 5. 
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TABLE 1.- BLADE CHARACTERISTICS 


Rotor radius, R 

8.53 m 

Blade lift slope, a 

5.73 

Blade semi-chord, b 

0.2185 m 

Number of blades, N 

4 

Rotor rotational speed, 0, 

23.67 rad/sec 

Blade twist rate, 0^ 

-8° 

Position of flapping hinge, rg 

0.3 m 

Blade cutoff, vq 

0.594 m 

C.G. position of blade, r^c 

2.74 m 

Position of lag damper, r^ 

0.3 m 

Lag damper coefficient, C ^ 

1000.0 N-m-sec/rad 

Blade mass, m 

106.4 kg 

Moment of inertia of blade, Ig 

1593.9 kg-m 2 

Mass moment of blade. Mg 

1659.1 kg-m 2 

Lock number, y 

8.84 

Wing section 

NACA 0012 

Gross weight, W 

62259.4 N 


TABLE 2.- BLADE PROPERTIES 


I 

Length, 

m 


EI 71 

N-m 2 

xlO 4 

EI Z , 

N-m 2 

xlO 5 

AI, 

kg-m 2 

xlO -3 

Am, 

kg 

xlO -1 

A6, 

deg 

xlO -1 

e, / 
m 

xlO -2 

1 

0.003 

34 

13 

22. 

37 

26.68 

94.71 

3.063 

0.0 

0.0 

2 

.327 

34 

13 

143 

4 

13.19 

87 ..80 

320 . 7 

.0 

.0 

3 

.254 

12 

97 

20 

28 

6 . 

801 

59.35 

144.9 


. 

0 

4 

.508 

6 

017 

4. 


4. 

929 

29.81 

24.06 

-.28 


4167 

5 

.610 

5. 

018 

4. 

214 

4. 

567 

31.57 

38.50 

-6.81 



6 

.280 







40.37 

24.21 

-3.12 



7 

.787 







47.02 

62.13 




8 

.457 







47.82 

39.67 

mSM 



9 

.737 







43.22 

58.05 

-8.23 



10 

.406 







44.71 

35.29 

-4.54 



11 

.711 







55.68 

53.82 

-7.94 



12 

.762 







38.35 

57.75 

-8.52 



13 

.254 







34.55 

22.02 

-2.83 



14 

.610 







47.02 

48.13 

-6.81 



15 

.635 







39.70 

60.23 

-7.09 



16 

.356 







23.03 

29.17 

-3.97 



17 

.254 





J 


33.22 

30.04 

-2.84 

J 


18 

.280 


r 

1 

i 


r 

7.448 

22.02 

-3.12 
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TABLE 3.- ROTOR OPERATING CONDITIONS 


Advance ratio, y 

0.18 

Collective pitch angle, 0 O 

8.0° 

Longitudinal pitch angle, 0 lg 

-1.48° 

Lateral pitch angle, 0 lc 

.695° 

Inflow ratio, A 

.0179 


TABLE 4.- TRANSFER FUNCTION OF ROTOR FOR GUST FREQUENCY OF 2.0 Hz 




-25.52 

-576.5 

-192.6 

22.40 



576.5 

-25.52 

22.40 

-192.6 



9.831 

-59.62 

-211.6 

-44.92 



66.51 

15.42 

3.635 

-179.3 



34.46 

-112.7 

2397. 

-85.89 



112.7 

34.46 

85.89 

2397. 

[T C ] 

= 

-63.17 

2.580 

-48.63 

216.3 



-2.580 

-63.17 

-216.3 

-48.63 



-1185. 

4.166 

96.97 

2508. 



-4.166 

-1185. 

-2508. 

96.97 



969.5 

77.35 

-10.06 

-293.7 



-77.35 

969.5 

293.7 

-10.06 

Advance 

ratio = 

= 0.18 



Rotor 

rotational speed 

= 23.67 rad/sec 

Gust 

frequency 

= 2.0 Hz 


Gust amplitude: 

U Go 

V Go ‘ °-° 

’ W Go 

{z} = 

[T c ]{0> 




{z} = 


: C’ T S’ 

He* H S’ 

M YC’ M YS’ 

Y C . Y 

{0} - 

(S 

>1’ 9 2’ 

6 3’ V 

V V T 



81.50 

-11.92“ 

11.92 

81.50 

50.86 

53.35 

-88.02 

59.16 

2561. 

91.49 

-91.49 

2561. 

-56.87 

37.14 

-37.15 

-56.87 

87.10 

-2618. 

2618. 

87.10 

4.302 

202,5 

-202.6 

4.305_ 


= 1.0 m/sec 
S’ M XC’ M XS’ M ZC» 







TABLE 5.- TRANSFER FUNCTION OF ROTOR FOR STEP GUST 


[T C ] = 


" 0.996 8.545 2.984" 

-0.496 8.004 2.077 

1.849 1.872 0.635 

.68.42 252.3 95.31 

-16.62 289.0 81.92 

. 8.599 17.79 5.652. 


Advance ratio =0.18 

Rotor rotational speed = 23.67 rad/sec 
Gust frequency = 0.0 Hz (step gust) 

Gust amplitude: u^ = Vg o = 0.0, w Gq = 1.0 m/sec 

{z} = [T c ]{©> 

{z} = (AT 0 , AH 0 , AMy 0 , AY q , AM Xo , AM Zq ) T 

{0} = ( 0 x » 0 2 , 0 3 ) T 
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COORDINATE SYSTEMS 


ROTOR: X R ,Y R ,Z R 
HUB: X H ,Y H ,Z H 

BLADE: x , y , z 



SINUSOIDAL GUST 


REFERENCE 
BLADE 
Xu 


(a) Gust shape and coordinate systems. 


Z R> Z H 



^ x 

X R 


(b) Blade deflection geometry. 

Figure 1.- Coordinate systems and blade deflection. 
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Figure 3.- Hub forces and moments. 
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(a) Horizontal force, F x . 


Figure 4.- Rotor impedance to gust excitation without feedback: w G = 1,0 m/sec 
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(b) Side force, F . 
Figure 4.- Continued. 
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NONDIMENSIONAL GUST FREQUENCY 


(d) Rolling moment, 
Figure 4.- Continued. 






PHASE ANGLE, deg AMPLITUDE, N/m/sec 




ojq/S1 

NONDIMENSION AL GUST FREQUENCY 

(e) Pitching moment. My.. 
Figure 4.- Continued. 
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NONDIMENSION AL GUST FREQUENCY 

(f) Yawing moment, M z . 
Figure 4.- Concluded. 






REAL 

(a) A0 = • Aw. 

Figure 5.- Stability analysis (root loci): y = 0.18 
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REAL 


(c) A0 = • Aw. 

Figure 5.- Concluded. 
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AMPLITUDE OF OSCILLATORY VERTICAL FOF 



0 -.2 -.4 -.6 -.8 - 1.0 - 1.2 - 1.4 - 1.6 

FEEDBACK GAIN, K w , rad/m 


(a) Aw feedback, (j> = 0°, 

Figure 7.- Sensitivity of the feedback gain to the vertical hub force, sensor at 
blade tip position: w~ = 1.8 m/sec, f r = 0.5 Hz, y = 0.18. 


AMPLITUDE OF OSCILLATORY VERTICAL FORCE 



FEEDBACK GAIN, K v X 10 2 , rad/m 

(b) Av feedback, <j> = 30°. 

Figure 7.- Continued. 


AMPLITUDE OF OSCILLATORY VERTICAL FORCE, 


X 10" 4 X 10 3 



(c) A<j> feedback, cj> = 270°. 
Figure 7.- Concluded. 
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AMPLITUDE OF OSCILLATORY HUB FORCES 



AMPLITUDE OF OSCILLATORY HUB FORCES 



AMPLITUDE OF OSCILLATORY HUB MOMENTS 



AMPLITUDE OF OSCILLATORY HUB FORCES 



AMPLITUDE OF OSCILLATORY HUB MOMENTS 



Figure 9.- Time history 
sinusoidal gust: w 
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iponses without control for a vertical 
: G = 0.5 Hz,' y = 0.18, 0 = 226 rpm. 
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Figure 10.- Time history of hub gust responses with individual blade pitch control for 
a vertical sinusoidal gust: w Go =1.8 m/sec, f G = 0.5 Hz, y = 0.18, fl = 226 rpm, 

= -0.5 rad/m. 
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Figure 12.- Time history of hub gust responses with individual blade pitch control for 
a vertical sinusoidal gust: Wq 0 = 1.8 m/sec, fg = 0.5 Hz, p = 0.18, fi = 226 rpm, 

^ = -0.5 rad/m. 
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ON-LINE 

Figure 15.- Four regulator options. 



SELF-TUNING REGULATOR 


(a) Open-loop regulator. 

Figure 16.- Adaptive regulators for gust alleviation. 
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SELF-TUNING REGULATOR 


(b) Closed-loop regulator. 
Figure 16.- Concluded. 



Figure 17.- Sequence of controller implementation. 1 
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Figure 18.- Time history of hub gust responses without c 
sinusoidal gust: Uq q = 3.0 m/sec, fg = 2.0 Hz, y = 


•1 


ol for 
, Q = 2 
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Figure 19.- Time history of hub gust responses without control for a lateral 
sinusoidal gust: v Go =2.0 m/sec, f Q = 2.0 Hz, p = 0.18, fi = 226 rpm. 
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Figure 20.- Time history of hub gust responses without control for three dimensional 
gust: uq 0 = 3.0 m/sec, vq q = 2.0 m/sec, Wq d = 1.8 m/sec, fg = 2.0 Hz, y = 0.18, 

0, = 226 rpm. 
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Figure 21.- Time history of hub gust responses with adaptive blade pitch control for 
a vertical sinusoidal gust. Adaptive open-loop regulator without caution (two 
inputs, two outputs): Wq 0 = 1.8 m/sec, fg = 2.0 Hz, y = 0.18, fi = 226 rpm. 


83 








d) SIDE FORCE 




TIME, sec 


Figure 22.- Time history of hub gust responses with adaptive blade pitch control for 
a vertical sinusoidal gust. Adaptive open-loop regulator without caution (six 
inputs, 12 outputs): w Go = 1.8 m/sec, f q = 2.0 Hz, p = 0.18, = 226 rpm. 
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ADAPTIVE OPEN-LOOP CONTROL 



Figure 23.- Time history of the mean square thrust perturbation: fg = 2.0 Hz, 

V - 0.18. 
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Figure 24.- Time history of hub gust responses with adaptive blade pitch control for 
a vertical step gust. Adaptive closed-loop regulator without caution (three 
inputs, three outputs): Wq 0 = 1.8 m/sec, y = 0.18, ft = 226 rpm. 
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Figure 25.- Time history of hub gust responses with adaptive blade pitch control for 
a vertical step gust. Adaptive closed-loop regulator without caution (three 
inputs, three outputs): wq 0 = 1.8 m/sec, y = 0.18, fi = 226 rpm. 
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Figure 26.- Time history of hub gust responses with adaptive blade pitch control for 
a vertical step gust. Adaptive open-loop regulator without caution (three inputs, 
six outputs): w Go = 1.8 m/sec, y = 0.18, fi = 226 rpm. 
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Figure 27.- Time history of quadratic performance function J: w q q = 1.8 m/sec. 
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Figure 28.- Time history of hub gust responses with adaptive blade pitch control for 
a vertical sinusoidal gust. Adaptive open-loop regulator with caution (six inputs, 
12 outputs): w Go = 1.8 m/sec, f G = 2.0 Hz, y = 0.18, Q = 226 rpm. 
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Figure 29.- Time history of hub gust responses with adaptive blade pitch control for 
three dimensional gust. Adaptive open-loop regulator without caution (six inputs, 
12 outputs): Uq q = 3.0 m/sec, vq 0 = 2.0 m/sec, wq d = 1.8 m/sec, f q = 2.0 Hz, 

p = 0.18, £2 = 226 rpm. 
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Figure 33.- Schematic view of lift and drag. 
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